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Abstract 

We investigate Ussing’s unidirectional fluxes and flux ratios of charged tracers 
motivated particularly by the insightful proposal of Hodgkin and Keynes on 
a relation between flux ratios and channel structure. Our study is based on 
analysis of quasi-one-dimensional Poisson-Nernst-Planck type models for ionic 
flows through membrane channels. This class of models includes the Poisson 
equation that determines the electrical potential from the charges present and 
is in that sense consistent. Ussing’s flux ratios generally depend on all physical 
parameters involved in ionic flows, particularly, on bulk conditions and channel 
structures. Certain setups of ion channel experiments result in flux ratios that 
are universal in the sense that their values depend on bulk conditions but not 
on channel structures; other setups lead to flux ratios that are specific in the 
sense that their values depend on channel structures too. Universal flux ratios 
could serve some purposes better than specific flux ratios in some circumstances 
and worse in other circumstances. We focus on two treatments of tracer flux 
measurements that serve as estimators of important properties of ion channels. 

The first estimator determines the flux of the main ion species from measurements 
of the flux of its tracer. Our analysis suggests a better experimental design so 
that the flux ratio of the tracer flux and the main ion flux is universal. The second 
treatment of tracer fluxes concerns ratios of fluxes and experimental setups that 
try to determine some properties of channel structure. We analyze the two widely 
used experimental designs of estimating flux ratios and show that the most widely 
used method depends on the spatial distribution of permanent charge so this flux 
ratio is specific and thus allows estimation of (some of) the properties of that 
permanent charge, even with ideal ionic solutions. The work presented in this 
paper is a first step showing how measurements of fluxes and flux ratios can give 
important insights into channel structure and function. 
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1 Introduction 


Challenges in study of ion channel functions: The atomic structure of the chan¬ 
nel protein is important for functions of ion channels (see, e.g., [9, 10, 16, 17, 18, 22, 27, 
30, 36, 37, 43, 44, 60, 61, 64, 71, 74, 80, 81, 82]). The atomic structure determines the 
class of mechanisms of ion channel functions, in an important but qualitative sense, 
and provides the basis for describing correlations in continuum models, increasing 
the resolution and realism of these models in significant ways. The structures being 
discussed are determined from crystals, with somewhat perturbed locations of atoms 
and forces (or they would not crystallize), in solutions remote from the physiological 
solutions in which the channels function, without gradients of electrical or chemical 
potential essential for channel function, and often at temperatures around 100K. The 
disorder of crystals is much less at 100K than at the biological temperatures (300K, 
more or less) where they were evolved to function. Entropic terms important at 
300K are much less visible at 100K. Changes in entropy with temperature can change 
the qualitative properties of materials. Semiconductors are insulators at low tem¬ 
peratures when their entropy is small but conductors at higher temperatures where 
entropy is larger. The structural approach does not deal at all with experimental 
measurements of current (that are the biological function of most channels). Current 
vs. voltage or concentration curves are not found in papers using the qualitative 
approach to selectivity. For example, a recent summary ([65]) in a leading journal 
includes neither current voltage relations, nor equations. Qualitative study is of great 
importance but it must not be confused with attempts at quantitative models. It is 
difficult to compare or improve qualitative models because they do not provide a 
basis for comparison. Qualitative models omit the specific fits to experimental data 
that allow one model to be compared to another. Quantitative models are needed to 
provide such comparisons. 

The quantitative treatment of ion channel functions and structures taken in this 
paper (reviewed in [11, 12, 13, 14, 19, 20, 23, 24, 25, 26, 31, 33, 34, 35, 47, 48, 53, 66, 
68, 69], etc.) differs from the usual classical treatment in several important respects. 
Most importantly the quantitative approach used here includes structure. Classi¬ 
cal treatments omit most structural information including the structure of electrical 
charge (permanent and polarization) that underlies so much of protein function and 
is canonized with the classical biochemical classification of amino acid side chains, 
as acid (negative permanent charge), basic (positive permanent charge), nonpolar 
(small dielectric coefficient) or polar (large dielectric coefficient). The quantitative 
approach based on continuum models treats ion channel functions as expressions of 
the current flow through a channel under a variety of conditions, namely different 
membrane potentials, different concentrations and compositions of ionic solutions, 
and different mutations. The quantities from computations/analyses of models can 
be compared directly with experimental measurements of current. The quantitative 
models are dramatically reduced in complexity compared to structures or simulations 
of structures in atomic detail, but they are precise. Such is the nature of most physical 
models of condensed phases. It is obviously impossible to include all the interactions 
of the tremendous number of particles involved. Both qualitative and quantitative 
approaches are needed and complement each other. 

It is convenient to treat ion channels as electric nano-devices with well defined 
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reasonably robust input output relations and easily predicted interactions with their 
surroundings. The quantitative treatments of ion channel functions are extremely 
important due to nanoscale, multi-scales, nonlocal interaction, limitations of experi¬ 
mental measurements, measurable characteristic quantities of ionic flows vs. charac¬ 
teristics of channel structures. 

Importances of a mathematical theory and numerics: The challenges in un¬ 
derstanding ion channel functions and structures strongly suggest the importance 
of mathematical analysis and numerical simulations as complementary tools to the 
physiological theory and experiments. A sound mathematical theory provides direct 
relations between the outputs and the inputs, at least for the simplified settings used 
in actual biological experiments. Once tested experimentally, the theory can make 
predictions and guide numerical designs for simulations and actual experiments. 

The enormous literature on channels and transporters from 1947 to say 1990 used 
measurements of fluxes, fluxes of radioactive isotopes, and their ratios as a main tool. 
The purpose of this paper is to define such fluxes in precise mathematical language so 
modern theories can be compared with those thousands of papers. Modern theories 
compute potentials from distributions of charge and have the advantage of being 
compatible with the Maxwell description of electrodynamics. Modern measurements 
are mostly of the structure and fluxes of channels and transporters. We believe 
measurements of structures, and fluxes, must be combined with electrodynamics to 
understand how channels and transporters work. We begin this process here by 
taking one precisely defined class of models (PNP type) and defining and computing 
the fluxes and flux ratios measured with such effort in so many laboratories for so 
many years. Transporters and channels have been studied in thousands of papers 
because they control most biological functions and are directly involved in a wide 
range of diseases. Our approach will identify the issues involved and pave the way for 
computations of fluxes and flux ratios in the substantial sequence of models that seek 
to remove the obvious limitations of classical PNP. It is the purpose of this paper to 
examine the Poisson-Nernst-Planck theory to see if the consistency of these models 
is significant. It turns out the consistency of the PNP type models is crucial. 

1.1 Important previous works on flux ratios 

The arguably most important quantity for ion channel functions is the total current, 
which is nowadays routinely measured based on an electric property of ionic flows. 
The total current is a combination of individual fluxes of all ion species in the steady 
state where, in particular, displacement current ( eodE/dt where E is the electric 
field) is zero, and provides crucial information on ion channel properties such as 
permeation and selectivity. Of course, the fluxes of individual ion species will provide 
more information. While individual fluxes are (additive) components of the total 
current, they cannot be decomposed from measurements of the total current using 
just electrical measurements. The flux ratio is needed as well, together with the flux 
measurements based on radioactive properties. 

Ussing’s work on flux ratios. Ussing, a physiologist, considered fluxes soon after 
the second world war, when radioactive isotopes of sodium first became available. At 
that time the role of proteins in membranes was not known so of course his analysis 
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could not include effects of the structure of proteins let alone their spatial distribution 
of charge. Indeed, the modern idea of protein structure was then just a dream in the 
mind of a few young workers mostly at Cambridge (UK and USA). In [79], based 
on an insightful observation from the Nernst-Planck equation (without using Poisson 
equation for electric potential explicitly), Ussing introduced and studied the flux ratio 
between the influx ,/M (the inward unidirectional flux from outside x = 0 to inside 
x = 1 of a channel measured at x = 1) and the efflux (the outward unidirectional 
flux measured at x = 0) of the same ion species. Note that the channel length is 
normalized to 1 (see (5.2) for a normalization). Ussing’s formula (equation (8) in 
[79]) is recast below 


J_lJM c M e ^o 

Jio] /(1) c(l) 
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c(l) 


( 1 . 1 ) 


where the quantity eo is the (positive) elementary charge, ks is the Boltzmann con¬ 
stant and T is the absolute temperature, z is the valence (number of charges per 
ion) of the ion species, c(0) and c(l) are the concentrations of the ion species at left 
boundary x = 0 and the right boundary x = 1, respectively, fi ex ( 0) and p ex (l ) are 

1 LL eX (0) 1 LL eX ( 1) 

the respective excess components (so /(0) = e k s T and /(1) = e k s T are 
the respective activity coefficients), Vo is the transmembrane electic potential - the 
difference between the electric potential at x = 0 and that at x = 1. 

An apparent advantage of flux ratios over individual fluxes is that its dependence 
on physical parameters involved in ionic flow is much simpler. As far as we know, 
physical scientists had not anticipated the importance and subtleness of the ratio of 
unidirectional fluxes despite the large physical literature on Nernst-Planck equations. 
Physiologists seem to have been the first to recognize the physical importance of 
flux ratios. This essentially abstract, physical and mathematical estimator was used 
universally to distinguish transporters from channels in thousands of experimental 
papers from 1949 to 1990 or so. The estimator proved to be quite robust in classical 
models that assumed electric potentials independent of the distribution of charge 


([5, 6, 62, 63]). 

Note that formula (1.1) provides a direct relation between the flux ratio with the 
data at two points x = 0 and x = 1. On the other hand, the activity coefficients 
/(0) and /(1) (or equivalently, the excess components p ex (0) and p ex (l)) cannot 
be determined by the data (Vo, c(0), c(l)) alone. Multicomponent systems like ionic 
solutions are complex fluids in which ‘everything depends on everything else’, as 
shown directly in the thousands of experimental measurements of activity coefficients 
in different solutions. The electrical potential is a global quantity depending on (for 
example) boundary conditions far away. The concentrations (more generally the 
chemical potentials) can depend on the shape of the structure confining them and on 
the electrical potential as well as on the interactions with ions of different type and 
variable concentration. 

Teorell ([76]), and Hodgkin and Huxley ([40]) used the ‘independence principle’ 
that states, with a constant potential difference across the membrane, the chance of 
any individual ion crossing the membrane in a given time interval is not affected by 
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the other ions that are present, to derive the flux ratio formula below 


— = C (°) c z CqVq 
JH c( 1) ’ 


( 1 . 2 ) 


where Co = The independence principle is useful, we now know, for highly 

selective channels that function independently in membranes (like the sodium and 
potassium channels studied by Hodgkin and Huxley). The independence principle 
does not apply to less selective channels (of which there are many, including the 
acetylcholine activated channel studied extensively by Katz, a collaborator of Hodgkin 
and Huxley) or to crowded bulk solutions. 1 

Formula (1.2) is simpler than (1.1). In particular, the right-hand side of formula 
(1.2) involves only the data (Vo, c(0), c(l)) at the boundary points. But formula (1.2) 
requires the independence hypothesis, and hence, could only be used to compare with 
experimental measurements of flux ratios as a TEST for independence hypotheses. 

Work of Hodgkin and Keynes. In [45], Hodgkin and Keynes pointed out that 
“Examples of mechanisms to which such an independence relation should apply are 
those involving combination with carries where only a small proportion of the car¬ 
rier is combined, or systems involving diffusion where the ions are present at fairly 
low concentrations throughout the membrane, and have no tendency to concentrate in 
narrow channels .” In comparing their experimentally measured flux ratio from potas¬ 
sium channels with that from formula (1.2) for a broad range of experimental designs, 
they found the deviation between these two flux ratios could be significant and con¬ 
cluded that independence hypotheses do not hold for ion channels in general. It is a 
fascinating fact that flux ratios for sodium channels (which were of great interest to 
Hodgkin and coworkers) have very different properties from potassium channels and 
do NOT show the interactions so characteristic of potassium channels. In fact sodium 
channels follow the independence principle ([7, 72]) nearly exactly, while potassium 
channels do not. Perhaps sodium channels function at very low occupancy and so 
only one ion is in the channel at a time and that ion cannot interact with neighbors 
even though the channel is narrow ([7, 15, 72]). 

In the same paper ([45]), based on a fitting function to the experimental data, 
Hodgkin and Keynes introduced the flux ratio exponent n' through the relation 
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(1.3) 


They argued that the flux ratio exponent n' estimates the average number of sites 
occupied by ions in the membrane. There is an enormous biological, physical, chemi¬ 
cal and mathematical literature building on the original work of Hodgkin and Keynes 

1 The independence principle played a large part in the early thinking of Hodgkin, Huxley, and Katz 
([38, 39, 40, 41, 42, 43, 44]) and was widely used by physiologists and biochemists to describe bulk 
solutions perhaps for that reason. Hodgkin was unaware of the evidence (personal communications 
ALH to RSE) that the Kohlraush principle of independent migration (~ 1880) had been disproven 
by measurements showing that properties of dilute sodium chloride solutions depended on the square 
root of concentration (because of the screening of the ionic atmosphere approximated by Debye- 
Hrickel theory). 
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on single file diffusion because nearly one dimensional systems are found in ion chan¬ 
nels, zeolites, and throughout the physical and biological world (see, e.g., [8, 52] and 
references therein). Hille summarizes the biological situation in Table 15.1 in [36]. 

In practice, one could determine the flux ratio exponent n' easily: by computing 
the value ^jye 2 ^ 0 ^ 0 from experimental designs of boundary conditions and experi¬ 
mentally measuring the flux ratio and then using the relation (1.3). This 

procedure of determining the flux ratio exponent n' has been used in thousands of 
experimental papers. 

The treatment of Hodgkin and Keynes, and its users following Hille ([36]) more 
than anyone else, involved electric fields that did not depend on the charge of the 
protein. Proteins are highly charged and the structure of their charges is an important 
determinant of their function: the location of acid and base groups and polar and 
nonpolar amino acids is a structural feature of the greatest importance. Charge 
changes electric fields in all models of electrodynamics (e.g., Maxwell equations) so 
the question arises whether the Hodgkin-Keynes-Hille interpretation of the flux ratio 
exponent is correct. The interpretation of n! by Hodgkin and Keynes in [45] was a 
startling insight and was widely thought to be a general property of channels resulting 
from their narrow structure (see Hille [36]) even though sodium channels in fact turn 
out to have independent fluxes similar to those in bulk solutions. The question is how 
much can we rely on the conclusions of this experimental work interpreted by using 
a theory assuming an electric field independent of protein structure and charge? 

To understand the flux ratio of Hodgkin-Keynes-Hille, we need a theory that in¬ 
cludes the channel structure, mainly its permanent charges. Analysis could determine 
how the crucial property of the flux ratio varies with conditions. The simplest treat¬ 
ment of channels in which the electric field is computed from charge is probably the 
classical Poisson-Nernst-Planck theory (see, e.g., [1, 2, 3, 4, 28, 29, 56, 55, 59, 70]). 
The classical PNP is just a beginning since it ignores correlations produced by the 
crowding of ions of finite size, although it does include long range correlations pro¬ 
duced by charge through the Poisson equation. We study flux ratios of the PNP type 
that also includes ion-to-ion correlations due to ion sizes, trying to present a formu¬ 
lation that can easily accommodate generalizations that include more correlations. 

Universal flux ratios and specific flux ratios. It seems that formula (1.1) sug¬ 
gests that flux ratios depend on bulk conditions only. It is not the case in general. 
First of all, Ussing’s design that leads to his formula (1.1) is special in the sense that 
the two unidirectional fluxes are produced in one experiment (essentially the same 
setup as in the two-isotope-setup discussed in Section 4.2). Still, in this case, the ex¬ 
cess potentials fi ex (0) and ^i ex (l) at x = 0 and x = 1 are involved in the formula. The 
excess potentials /a ex (0) and /i ex ( 1) generally also depend on the channel structure 
except when we use local models (see Hypothesis 2.2) for excess potentials. Other 
experiment designs will lead to flux ratios that depend on channel structures even 
when local models are used for excess potentials, in fact, even if the classical PNP 
model is used where the excess potentials are totally ignored. 

It is reasonable and important to classify types of flux ratios due to different 
experiment designs, relative to the channel structure, into two types: universal and 
specific. More precisely, if the flux ratio from a certain experimental design is depen¬ 
dent only on properties of bulk solutions (boundary conditions), particularly, not on 
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the channel structure, the flux ratio is called universal (see part (i) and Setup 2 in 
part (ii) of Section 1.2 for examples of universal flux ratios); if the flux ratio depends 
on the channel structure too, then it is specific (see Setup 1 in part (ii) of Section 1.2 
for an example of specific flux ratios). 

We emphasize that it is the detail of an experimental design that dictates the type 
of flux ratios. For the same channel, some experimental designs could lead to universal 
flux ratios (see Corollary 3.2 and particularly Remark 3.3); some experimental designs 
could lead to specific flux ratios (see Corollary 4.2). 

We point out that this classification is rough but helpful. The flux ratio in for¬ 
mula (1.1) is, strictly speaking, specific. But, for reasonably dilute bulk solutions, 
local model for excess potentials are good enough so the flux ratio is nearly universal. 
Also, the terminology could be misunderstood to mean that the universal flux ratios 
are useless. Indeed, universal flux ratios are not directly useful for detecting chan¬ 
nel structures. But universal flux ratios are helpful in understanding other channel 
properties (see the discussion of treatment (i) of flux ratios in Section 1.2). 

1.2 Experiments on flux ratios and Poisson-Nernst-Planck theory 

Motivated by the works of Ussing, and of Hodgkin and Keynes, flux ratios have been 
used in many studies of various properties of ion channels. Two major applications 
will be briefly described below and will be carefully analyzed based on PNP type 
models in the main part of this paper. 

(i) Flux ratios of a main ion species and its tracer. One application of flux ra¬ 
tios is to determine the flux of a main ion species by measuring the flux of its tracer 
(a radioactive isotope). Roughly speaking, one adds a small amount of a radioactive 
isotope of the main ion species and measures the flux of the isotope from its radioac¬ 
tive property , namely the number of disintegrations per unit time (detected by the 
radiation products of the disintegration such as gamma rays), and then estimates the 
flux of the main ion species by determining the flux ratio - proportionality constant 
- between the flux of the tracer and that of the main ion species. In view of Ussing’s 
work, the proportionality constant can be well approximated and easily determined 
from the boundary conditions alone (see Corollary 3.2). 

An ideal experimental design for this purpose uses a universal flux ratio so that 
the proportionality constant can be precisely and easily determined, independent of 
the often unknown channel structure. Based on analysis of PNP type models, we 
discover an experimental design (proposed in Remark 3.3) so that, if the design can 
be practically implemented, the resulting flux ratio is indeed universal. 

(ii) Flux ratios vs. channel structures. Another setup for flux ratios is moti¬ 
vated by the work of Hodgkin and Keynes ([45]) discussed above. These setups are 
designed so the estimated flux ratios are NOT universal. Rather these flux ratios sup¬ 
pose to depend on the characteristics of the channel, which would allow one to extract 
information on channel structure - permanent charge, in principle. We consider two 
different experimental setups (see Section 4) that produce estimates of flux ratios. 
One setup, the one-isotope-setup (Setup 1), uses one isotope in two experiments, and 
the other, the two-isotope-setup (Setup 2), uses two isotopes in the same experiment 
(essentially the same setup Ussing used in deriving formula (1.1)). With the specific 
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setups described in Section 4, it will be clear that the two flux ratios from the two 
setups are nearly, but not quite equal (numerically). 

Analysis of PNP type models indicates that the issue of experimental setups in 
this application of flux ratios is subtle. In Section 4, we analyze the flux ratios 
from the two different setups based on PNP type models with local excess potential 
models. Our analysis shows that the two different setups make a significant qualitative 
difference - the flux ratio from the one-isotope-setup (Setup 1) is specific so the flux 
ratio depends on the permanent charge, even with ideal ionic solutions; the flux ratio 
from the two-isotope-setup (Setup 2), for PNP type with local excess potential models, 
is universal, totally independent of the permanent charge. 

The ultimately important issue is the relation between the flux ratio and the 
channel structure. As mentioned above, for the two-isotope-setup (Setup 2), the 
PNP model with local excess potential does not provide any relation between the 
flux ratio and the channel structure. But, for the one-isotope-setup (Setup 1), PNP 
with local excess potentials, even PNP with only ideal electrochemical potentials , does 
provide such a relation although not in an explicit form yet. To further examine the 
relevance of the flux ratio from the one-isotope-setup (Setup 1) to channel structures, 
we consider the case where the permanent charge is small relative to the characteristic 
concentration and a local hard-sphere is included in the electrochemical potential. An 
approximation formula of the flux ratio is obtained with the leading terms explicitly 
expressed in terms of boundary conditions and the assumed channel structure. This 
formula allows one to analyze the relation between the flux ratio and the permanent 
charge in this simple case. We emphasize that our work in this paper is only the 
starting point of our investigation so the results are far from a complete theory on re¬ 
lationship between the flux ratio and the channel structure. Indeed, one of the reasons 
for this paper is to pose this important biological language in precise mathematical 
form so more complete theories can be developed and applied. 

A comment on the sign of unidirectional flux. In formulas (1.1), (1.2) and (1.3), 
./M and ./M are both defined as positive quantities, although they flow in opposite 
directions. In the following, we will fix one spatial orientation in our formulas for all 
fluxes. Thus jM is positive and j[° 1 is negative. This produces a minus sign in our 
formulas for flux ratio jW/jH. 

Organization of the paper. In Section 2, a quasi-one-dimensional PNP type model 
for ionic flow is reviewed and our assumption on the excess potentials is introduced. 

In Section 3, a formula (3.1) for individual fluxes is provided. We then examine the 
proportionality between the flux of the main ion species and the unidirectional flux 
of a tracer. A formula (3.7) for the proportionality constant is derived and, based on 
the formula, we propose an extremely useful experimental design in Remark 3.3 from 
which this flux ratio is universal - the ideal situation for the purpose of determining 
the proportionality constant. 

In Section 4, we discuss two widely used experimental setups for measuring unidi¬ 
rectional fluxes and their flux ratios of unidirectional fluxes, and their differences are 
discussed. For PNP type models with local excess potentials, the analysis indicates 
that the flux ratio from the two-isotope-setup (Setup 2) is universal and does not 
provide any information on channel structures but the flux ratio estimated from the 
one-isotope setup (Setup 1) is specific and contain information of channel structure. 



In Section 5, we consider a special case for the one-isotope-setup (Setup 1) to 
illustrate that the flux ratio is indeed specific , in a concrete way. The permanent 
charge is described in (5.1) and is assumed to be small relative to the characteristic 
concentration, a local hard-sphere approximation for excess potential is specified in 
(5.5), and, particularly, the dimensionless quantity e defined in (5.2) is assumed to 
be small. We only examine the zeroth order (in e) term of this specific flux ratio. 
(For convenience, we refer to this zeroth order (in e) term of the flux ratio simply as 
the flux ratio.) Then approximation formulas (5.14) and (5.18) for the flux ratios are 
derived for the one-isotope-setup (Setup 1) with different boundary electroneutrality 
conditions. The approximate formulas for the flux ratios are expressed explicitly 
in terms of boundary conditions and the prescribed permanent charge distribution. 
Section 6 contains a summary of our results and further discussions on our work. 

The derivation of formulas (5.14) and (5.18) for the one-isotope-setup (Setup 1) 
is provided in the Appendix (Section 7). Section 7.1 gives a reformulation of the 
flux ratio formula for the special case considered in Section 5. Formula (5.14) for 
electroneutrality boundary conditions among all three ion species is derived in two 
parts: Section 7.2.1 handles the contribution to the flux ratio from (small) permanent 
charge and Section 7.2.2 treats the contribution of ion sizes with a local hard-sphere 
potential. Formula (5.18) for electroneutrality boundary conditions among only the 
main ion species and the counter ion species is derived based on the derivation of 
formula (5.14) with necessary modifications indicated. 

2 PNP type models and electrochemical potentials 

2.1 The boundary value problem 

The following quasi-one-dimensional PNP type model was suggested by Nonner and 
Eisenberg [67] (see [57] for a justification from the three-dimensional PNP for a special 
case): 


^4( erWeoA(x) jt^) = + <?(*)), 

v ' e— I 

d 

~1 Jk — 0, J k — 
dx 

where 0 is the electric potential, ks is the Boltzmann constant, T is the absolute 
temperature, eo is the elementary charge, e r (x ) is the relative dielectric coefficient, £q 
is the dielectric constant of vacuum; h(x) is the cross-section area over the longitude 
coordinate x of the channel, Q(x) is the permanent charge; for the kth ion species, c k 
is the concentration, z k is the valence, J k is the flux density across each cross-section 
of the channel, D k (x) is the diffusion coefficient, fi k is the electrochemical potential 
(that depends on 0 and {cj} in general). 

The boundary conditions are, for k = 1,2,... ,n, 

<t>(0) = V o , c fc (0 ) = L k ; 0(1) = 0, c k (l) = R k . (2.2) 

The concentrations of ions manipulated in experiments are represented by the set of 
boundary values of concentrations L k s on the left side and R k s on the right side. 


——D k (x)h(x)c k —/i k , k = 1,2,... ,n 
ks-L dx 
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We remark that the channel has been normalized from x = 0 to x = 1 by a scaling of 
its longitude variable (see (5.2) for the full dimensionless rescaling). 


2.2 Electrochemical potential 

The electrochemical potential //*, = /if 1 + fif x , consisting of the ideal component 
and the excess component is the most important variable in a model of ions in 
solutions or channels. The ideal component 

= k B T In — + z k e 0 (t>, (2.3) 

co 

where Co is a characteristic concentration, reflects the contribution of ions as point- 
charges. For nearly infinite dilute ion solutions, the ideal component /rl d is a good 
approximation of the electrochemical potential /i/ c . 

The excess component f/if measures the deviation of the electrochemical potential 
from the ideal component, in particular, it accounts for ion size effects and presum¬ 
ably channel size in turn. This description of nonideal properties uses the excess 
component as an additive component to the ideal component. Another equiv¬ 
alent description of nonideal properties is in terms of activity coefficient (discussed 
below) as multiplication factor to the one associated with the ideal component. The 
nonideal property of kth ion species depends, in principle, on everything involved 
such as the temperature, the medium (water), ion-ion interactions due to ion sizes, 
the concentrations of all ion species of the mixture, etc. Models of are not fully 
understood and developed. There are various explicit approximations of /j| x that 
are particularly important in models of ionic solutions in and around biological cells. 
These are all derived from seawater, and so are quite concentrated mixtures (~ 0.2 M) 
and so have substantial and important excess components. Indeed, animal cells de¬ 
pend on the difference between potassium and sodium: without that difference they 
cannot survive, and potassium and sodium ions differ only because of their excess 
properties due to their different sizes. 

The activity a k is another standard characteristic of the electrochemical potential 
of kth ion species and describes nonideal properties through its activity coefficient f k 
as a multiplicative factor to the one associated with the ideal component (up to a 
reference potential). The activity a k is dimensionless and is defined by 


a k 


e k B rh l k Mfc) _ e k B T^k X e k B T^k Mfc) 


(2.4) 


where, in general, /i° can be any reference potential and, in this sense, the activity a k 
is a relative term. The activity coefficient f k is the pre-factor in the last expression 
in (2.4) 


/ fc = e *B^r (2.5) 

Note that f k = 1 for ideal ionic solutions. The deviation of f k from 1 measures the 
nonideal nature of the ionic solution. The relation (2.5) is used to get the second 
expression of the flux ratio in (1.1). For a solute in solution, a standard convention 
is to take = njf. In this case, one has a k = fk- If we take ^ = Zkeo4> (only the 
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electric component of ;T k ) and the characteristic concentration co = ffj Cj , then 

— fk “ — fh-I'h i 
CO 

where x k is the concentration fraction of fcth ion species, is another form of activity. 
If we take fi k = 0, then (2.4) can be rewritten as 

_1 C ]b 

a k = fk e ksT k or In a k = In f k + In — + z k ( 0 f. 

co 

The latter is the equation (3) with co = 1 in [79] used by Ussing to define the 
(electrochemical) activity a k that includes the electrical potential f> as a part. 

In the following, we will use the excess components yffs to characterize the 
nonideal nature of ionic solutions. 

For all results in Section 3 and Section 4, we will make two assumptions about 
the electrochemical potentials adopted in the PNP type model (2.1). 

Hypothesis 2.1. For a given fixed ionic mixture, if jth ion species and kth ion 
species have the same valence and diameter, then the excess potentials and diffusion 
coefficients for jth ion species and for kth ion species are the same; that is, if zj = z k 
and dj = d k , then Hj X (x) = pL e k (x) and Dj(x ) = D k (x). 

The validity of Hypothesis 2.1 in great generality was explicitly pointed out by 
Ussing in the last paragraph on page 45 of [79]. In our work, Hypothesis 2.1 will be 
applied to a main ion species and its isotopic tracer that have the same valence and 
diameter. 

For a solution (4>(x), C(x), J) of the boundary value problem (BVP) (2.1) and 
(2.2), we will thus use the following notation for gf k (x) to indicate explicitly its 
dependence on ( z k ,d k ): 


dT(x) = H e k x (x;z k ,d k ). 


( 2 . 6 ) 


Hypothesis 2.2. The excess potentials are taken to he local in the sense that, for 
kth ion species, its excess potential n k x (x; z kl d k ) at location x depends pointwise on 
concentrations C(x) = { Cj(x ) : j = 1,... ,n} of all ion species of the ionic mixture 
but only at the location x; that is, pL e ff(x;z k ,d k ) = p e ff(C(x);z k ,d k ). 

For a nonlocal model of the excess potential pfff , its value pi e k {x) at x is often taken 
to be dependent on the values of C(y ) for y in a neighborhood of x where the size of 
the neighborhood is typically about the maximum of (dj + d k )/2 for j = 1,2,... ,n 
and the Debye screening length at x. 

The local nature of excess potentials assumed in Hypothesis 2.2 is an approxi¬ 
mation and is often good enough. We will further comment on this hypothesis in 
Remark 2.3 and at other places that it is crucial. 

In the following, we separate the ideal gas component from the electrical and 
nonideal terms to make our formulas more tidy and concise, that is, we introduce 


p k (x;z k ,d k ) 


Zk(o4>(x) + 


(x’j Z k ,d k ), 


(2.7) 
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where Co = eo/ksT. Therefore, for kth ion species, k B Tp k (x-, z k , d k ) is the component 
beyond the ideal gas (non-charged) component k B T\n{c k (x) / cq), that is, the total 
electrochemical potential p k (x) is 

/ifc(x) = 4 bT1ii^-^ + k B Tp k (x;z k ,dk). 
co 

We also purposely include z k and d k in the expression p k (x\ z k ,d k ) to remind the 
readers of Hypothesis 2.1; that is, for jth ion species and fcth ion species, if Zj = z k 
and dj = d k , then p k (x ; z k , d k ) = Pj{x ; Zj, dj ). 

Remark 2.3. We claim that, for all the results in Section 3 and Section 4, the local 
assumption (2.2) on the excess potentials p e k (x\ z k , d k ) ’s can he relaxed under one 
general condition; that is, one needs p1f(x-, z k , d k ) to he local only outside the channel 
(near x = 0 and x = l), in particular, nonlocal models for the excess potentials are 
allowed inside the channel, provided the corresponding BVP (2.1) and (2.2) is well- 
posed. This claim follows from the proofs of all the results in Section 3 and Section 
4 without any change. On the other hand, when nonlocal models for excess potentials 
are used at inside and outside (particularly near x = 0 and x = 1) of the channel, 
the corresponding BVP (2.1) and (2.2) is known to he severely under determined, 
and it is not completely understood how to formulate “correct” boundary conditions 
for the BVP to be well-posed. In [75], an initial attempt for a correct formulation of 
boundary conditions that are natural to the biological setups is taken. Further studies 
would be necessary for a better understanding of the roles of nonlocal nature of the 
excess potentials. 


3 Unidirectional flux and proportionality constant 

For the results in this section and in Section 4, we only assume that the BVP speci¬ 
fied in the equation (2.1) and the boundary condition (2.2) has a solution. It is not 
assumed that the solution of the BVP (2.1) and (2.2) is unique (of course, the results 
are referred to a given solution of the BVP (2.1) and (2.2)). We also do not assume 
electroneutrality boundary conditions so there might be boundary layers. Experi¬ 
mental designs often use “four electrode” arrangements to remedy possible boundary 
layers. In those arrangements, two ‘working’ electrodes (that carry current) are in¬ 
serted far away from the channel so boundary conditions that may produce boundary 
layers do not reach the ends of the channel. Two other electrodes are inserted near 
the ends of channel and away from possible boundary layers and used to measure 
electric potential and concentrations. 


3.1 An important formula for J k and unidirectional flux 


Proposition 3.1 below has been derived in classical channel models for ideal solutions 
where p k (x;z kl d k ) = z k ( 0 (f(x) (see, e.g., [49]). 


Proposition 3.1. Assume Hypotheses 2.1 and 2.2. If (f),C,J) is a solution of the 
BVP (2.1) and (2.2), then 


J k = ^L k e Pk{0 ’ Zk4k) - R k e Pk{1 ’ Zk ’ dk ^ 


pVk^i^k^k) \ —1 

__ dx I 

D k (x)h(x) ) 


(3.1) 
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where p k (x; z k , d k ) is defined in (2.7). 

In particular, if, at the boundaries (outside the channel), we do not include the 
excess component, this is, if we take p k ( 0) = iflfiO) an d hk{ 1) = 1 ^( 1 ), then 


J k = (L k e z ^ oVo - R k 



1 e Pk(x;z k ,d k ) N _i 

-dx 


/0 D k (x)hfx ) 

Proof. It follows from Nernst-Planck equations in (2.1) that 

- J k =D k (x)h{x)e~ pk ^ Zk4k) (c k e p ^ x ’ Zk ’ dk A'. 


(3.2) 


(3.3) 


Divide the three factors on the right-hand side of (3.3) and integrate to get 


Jk 


r l e Pk{x-,z k ,d k ) 

' 0 D k {x)h(x) 


dx = £ (c k e Pk{ - x ' Zk ’ dk ^dx = c fc ( l)e Pk ^ Zk4k) - c k (0)e Pki °’ Zk4k) . 


The formula (3.1) for J k then follows from the boundary conditions (2.2). 

When taking p k (0) = p 2 k (0) and p k ( 1) = p td (l), one has p k (0;z k ,d k ) = z k Co Vo 
and p k (l] z k , d k ) = 0. Formula (3.2) then follows from formula (3.1). □ 


Formula (3.3) is the same as Ussing’s formula (5) in [79], from which Ussing 
derived his formula (1.1). The extra step leading to formula (3.1) proves to be useful 
at several occasions, for example, for obtaining formulas (3.7) and (4.3) below which 
are not direct consequences of Ussing’s formula (1.1). 

Formula (3.1) for J k has the advantage that the only unknown quantity on the 
right-hand side is the function p k (x\ z k , d k ) which, from Hypothesis 2.1, is the SAME 
for ion species with the SAME ( z k ,d k ) such as for an ion species and its isotopes. 

We end this part by illustrating an immediate use of formula (3.1): a definition 
of unidirectional fluxes. This has been demonstrated for ideal ionic solutions in [49]. 

Formula (3.1) can be rewritten as follows to give a decomposition of J k as the 
sum of two unidirectional fluxes J k = j\} + where 


jf =L k e Pk{0 ’ Zk ’ dk) ( 

4° ] =-i4e Pfc{1; * fc ’ dfc )(^ 


1 e Pk( x ;z k ,d k ) N _i 


D k (x)h(x 


-dx 


p Pk{x-,Zk,d k ) . _1 

dx 


0 D k (x)h(x) 


(3.4) 


are the so-called influx and efflux of kth ion species, respectively. (Note the minus 
sign in jjf 1 is due to the fact that we use the same spatial orientation of positive 
current for both fluxes.) 

Formulas in (3.4) provide a mathematically meaningful definition of influx and 
efflux. It should be stressed that the unidirectional influx (respectively, efflux 

J [ ° ] ) depends also on all other quantities in the boundary conditions as well as the 
permanent charge of the problem through the global profile of p k (x; z k , d k ). The ratio 
between these two unidirectional fluxes 



J Jke Pk{ 0 -,z k ,d k ) 

R ke p k (l;zk,dk) 
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agrees with Ussing’s formula (1.1). When local excess potentials are used, the above 
ratio is independent of channel structures, and hence, the flux ratio is ‘approximately’ 
universal. This formula, under the hypothesis that, near x = 0 and x = 1, ionic 
solutions are dilute enough, would reduce to formula (1.2); in particular, the flux 
ratio is universal, independent of channel structures. 

3.2 Proportionality between flux of main ions and that of its tracer 

While most experiments today measure the total current, it is important to mea¬ 
sure the fluxes of individual ions - the components of the current - and such flux 
measurements were the basis of hundreds if not thousands of classical papers and 
almost all work on vesicle preparations widely used by molecular biologists. (The 
references [36, 37, 77] are good entry points to the classical literature.) For exam¬ 
ple, for Na + Cl - solutions, one can add a tiny amount of a radioactive isotope of 
sodium denoted by T^ a . The flux of T/v a can be measured by its radioactivity. The 
isotope T/v a has exactly the same electrical and chemical properties as the main ion 
species (the NON-radioactive Na) so one expects their fluxes are proportional to each 
other in a simple manner. Determining the proportionality constant would provide 
a measurement of the flux of the main ion species by measuring the tracer flux of 
radioactive T^ a . A general analysis of the proportionality constant or ‘specific ac¬ 
tivity’ (the favorite name of experimentalists) is not possible, of course, because the 
‘constant’ varies with experimental setup. Experimental setups leading to universal 
flux ratios are preferred. 

Since the value of the proportionality constant is different in different setups, we 
will consider one experimental setup that is widely used in the experimental literature 
and propose, in Remark 3.3, a specific experimental setup that has a substantial 
advantage. 

For simplicity, we consider three ion species in this subsection: the main ion 
species (for example, Na + ), its tracer (an isotope of Na + ), and the counter ion (for 
example, Cl - ). For the relative quantities, we use the following notation. We use 
the same subscript 1 for both the main ion species and its isotope to indicate that 
they have the same valence, etc., but use different superscripts [ m ] for the main ion 
species and [t] for the tracer to distinguish them. The symbols for the boundary 
concentrations are used in the similar way. 



Concentration 

Valence 

Flux 

Main ion species 

c i 

z\ > 0 

tH 

■ J \ 

Tracer 

c l<l 

C 1 

z\ > 0 

jV\ 

J 1 

Counter ion 

C2 

Z2 < 0 

■h 


For the moment, we do not specify the boundary concentrations of the tracer but 
will do so afterwards. Thus, the BVP (2.1) and (2.2) for the main ion species, the 
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tracer and the counter ion species becomes 


1 d 


h(x) dx v £ r(z)£°/i(a:)—0) = -e 0 ^i(c^ m] +c^) + z 2 c 2 + Q(x)^. 

d [m] _ d [t ] _d 

- J i ~~h 


dx 


- J mj = 


dx 

1 


d 


_ T W - 

1 fc B T 


— Di(z)h(x)<% ] 


“ * = T^f D2{x)h{x > c ^ 

with boundary conditions, 

m = Vo, 4 ml (0)=4 m] , cf ] (0) = Lf 1 , c 2 (0) = L 2 ; 
^(1) = 0, c H(i) = jR H, c W( 1 ) = j rM c 2 (1) = i? 2 . 


(3.5) 


(3.6) 


Recall that the main ion species and its tracer have the SAME valence, diam¬ 
eter and diffusion coefficient. The following result is then a direct consequence of 
Proposition 3.1, not a direct consequence of Ussing’s formula (1.1). 

Corollary 3.2. Assume Hypotheses 2.1 and 2.2. For BVP (3.5) and (3.6), the 

[tl fml 

proportionality constant between the flux J{ of the tracer and the flux J{ of the 
main ion species is 


jf ]jP] e pi(0;zi,<2i) _ _^W e Pi(l;zi,<ii) 

jj m l ^Hgpi(°;zi,di) _ ;z 1 ,d 1 ) ' 

Furthermore, if we take /xi(0) = pdfifl) and /xi(1) = 1), then 

jf ] _ L [ l ] e z ^ oVo - Rf 
“ L H e , lCo Vo _ i?H ‘ 


(3.7) 


(3.8) 


Proof. Since both the main ion species and the tracer have the same valence z\ and 
diameter d\ , and hence, the function pi(x; z\,d\) in (3.1) is the same for both of them. 
Also both the main ion species and the tracer have the same diffusion coefficient 
D i(x). Formula in Proposition 3.1 yields the formula for the quotient. □ 


Note that, since excess potentials are approximated by local models, the propor¬ 
tionality constant or the flux ratio given by formula (3.7) is approximately universal. 

We now apply the results to specific boundary concentrations of tracers widely 
used in experiments. It is arranged so that the radioactivity on one side of the channel 
is kept nearly zero (by flowing solution by the channel). If one takes > 0 and 
= 0, then (3.7) gives the proportionality constant 

jW _j,M e Pi(0;zi,di) 

j[m] xMgpROjzi.di) _ ^HgpRl^di) ’ 
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where ./[^ = jf’^ is the inward unidirectional flux (influx) of the tracer. This flux 
ratio, with any local models for excess potentials in the solution outside the channel, 
is universal. Given boundary conditions, one still needs to estimate p\(f)-, z\, d\) and 
pi(l; Z\, di), which depend on specifics of local models for excess potentials, to deter¬ 
mine the proportionality constant. If this experiment involves dilute ionic mixtures, 
then one can apply formula (3.8) to approximate the proportionality constant 

jf ] _ 

“ L M e * lf oVo - 4 m] ' 

Finally in the remark below, we make an observation and propose a possible design 
of experiments for a simple and precise determination of a proportionality constant. 

Remark 3.3. One experimental setup provides simple and exact results and so we 
recommend using it, where practicable. If one designs experiments so that 

Lf/L [ ™ ] = Rf ] /R [ ™ ] = a, (3.9) 


then it follows from formula (3.7) that 


■J 


[t] 




= a 


(3.10) 


where J = jf is the total flux of the tracer and j{”^ is the flux of the main 

ion species. Note that this simple relation is universal and remains true independent 
of specifics of local models for excess electrochemical components and even if nonlocal 
models are used as long as the boundary value problem is well-posed. 

If one can design this experimentally and measure the total tracer flux .l\^ = 
J-f’*] + jf '°\ then formula (3.10) gives us a beautiful result. It says that if the “specific 
activity” defined in equation (3.9) is the same in the left and right hand solutions, 
then the radioactive species does in fact trace the movement of nonradioactive species 
robustly essentially without error. We do not know whether or not this design has 
been discovered and used for tracer measurements and regretfully we do not know how 
to find the appropriate reference in the vast literature. 

There is a further powerful and beautiful result of this simpler and accurate for¬ 
mula. Consider a setup with two tracers T\ and Ttwo isotopes of the main ion 
species with different radioactivity. One designs the boundary concentrations for trac¬ 
ers as follows: 


for tracer T\ : >0 at x = 0 and R^ =0 at x = 1 ; 

for tracer T 2 : L ^ =0 at x = 0 and R ^ >0 at x = 1 . 

If L^/L^ = a and R^/R^ = a, then 


= - J 


a 


(jf 1 ’'] + 


where is the influx of tracer T\ and is the efflux of tracer T- 2 - Under these 

conditions, the tracer in fact traces the main ion flux exactly. 
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4 Flux Ratios from two setups 


In this section, we move to a different but related topic; that is, we examine the flux 
ratios of two unidirectional fluxes of tracers in two setups corresponding to two differ¬ 
ent experimental designs for the purpose of detecting channel structures. Two issues 
are relevant: the first (the focus of this section) is whether or not each flux ratio is 
universal or specific and is able to provide any information on channel structures, and 
to compare their differences based on PNP type models with local excess potentials; 
the second is the detailed relation between the specific flux ratios and the channel 
structure. The latter will be examined by case studies in Section 5. 

Since our PNP type models contain local excess potentials, our results apply for 
non-ideal ionic solutions and, as direct consequences, we also include relevant results 
for ideal ionic solutions. 


4.1 Flux ratio for the one-isotope-setup (Setup 1) 


For this setup, one uses one main ion species, one tracer (an isotope of the main ion 
species), and one counter ion species. Two sets of experiments are performed. 

In one experiment, the boundary conditions for the tracer are nonzero concentra¬ 


tion L^" 1 ' > 0 on left boundary x = 0 and zero concentration Rf’' 1 = 0 on the right 
x = 1. This would produce an influx of the tracer. 

In the other, the boundary conditions for the tracer are zero concentration = 
0 on left boundary x = 0 and nonzero concentration R^'°^ > 0 on the right x = 1. 
This would produce an efflux jf’°^ of the tracer. 

We will examine the flux ratio jf’^/ jf’°^ associated with these two experiments. 
According to the two different experiments, two different sets of boundary condi¬ 
tions are as follows. 


M 


(BVi) : 


# ] (0) = Vo, 
c l ™’ i] ( 0) = L [ ™’ i] > 0, 
4 M1 (0) = L [ l A > 0, 


0W(1) = 0, 

C M( 1 ) = R [rn,i] > 0j 

cf%) = r\ A = 0, 


c 2 (0) = 4 > 0, c 2 (1) = -^2 > 0- 


(main ions) 

(tracer) 

(counterions) 


(4.1) 


(BVo) : 


</> [ol (0) = Vo, 
c' m ’° ] (0) = 4 m ’ o] > 0, 
cf’ o] (0) = l\ A = 0, 

4 o] (0) = 4 o] > 0, 


0 [o] ( 1 ) = 0 , 

4 m ’ o] (i) = 4 m ’ o] > 0 , 


4’ OJ (i) = R l !’° l > 0, 
4 o] (i) = 4 o] > 0 . 


M 


(main ions) 

(tracer) 

(counterions) 


(4.2) 


We will thus consider TWO boundary value problems: one has system (2.1) with 
the boundary conditions (BVi) in (4.1) and the other has system (2.1) with the 
boundary conditions (BVo) in (4.2). The key difference between these two boundary 
conditions concerns the tracer concentration: it is zero at x = 1 for (BVi) and it 
is zero for (BVo) at x = 0. For simplicity, the boundary condition for the electric 
potential 4> will be kept the same for both BVPs. The boundary concentrations for 
the main ion species and those for the counter ion species could be different for the 
two BVPs for the purpose of detailed case studies in Section 5. 
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For the boundary condition (BVi) where L > 0 and = 0, the absorbing 

boundary is at the right end x = 1 where the tracer concentration is zero and the 

\t i] 

tracer is inserted from the left end with concentration L\ ’ 1 > 0 at x = 0; in particular, 

[t i] 

the flux of the tracer is the influx ’ . The notation in the next table will be used. 


BVi 

Concentration 

Valence 

Flux 

Main ion 

C 1 

zi > 0 

T [m,i\ 

Tracer 

c i 

z\ > 0 

j\ff\ 

Counter ion 

c i*j 

c 2 

22 < 0 

t!*J 

J 2 


Similarly, for the boundary condition (BVo) where L^'°^ = 0 and R^’°^ > 0, the 
flux of the tracer is the efflux jf'°^. For this case, we use the notation in next table. 


BVo 

Concentration 

Valence 

Flux 

Main ion species 

J«u°J 

c i 

z\ > 0 

j[m,o\ 

J 1 

Tracer 

Jfl°\ 

c i 

z\ > 0 

t\Vo\ 

J 1 

Counter ion 

c [ ° l 

c 2 

2 2 < 0 

J 2 


Concerning the flux ratio between j\ Ll ^ and jf’°\ we have the following formula. 

Theorem 4.1. Assume Hypotheses 2.1 and 2.2. For the one-isotope-setup (Setup 
1) with one tracer in two experiments, let be the influx of the tracer associated 
with the boundary conditions (BVi) in (4-1) for one experiment; let jf’°^ be the efflux 
of the tracer associated with the boundary conditions (BVo) in (4-2) for the other 
experiment. 

Then the flux ratio between and jf’is 


J 


M ]4tA(0\zi,<h) rl e p^(x; 2 i,di) 


dx 


J 


M R [t,o] e p [ fl(i-,z 1 ,d 1 ) Jo D 1 (x)h(x) \J 0 Di(x)h(x) 


"1 p [ f(x\zi,<h) 

- dx ) , (4.3) 


where py (x; z i, d\) defined in (2.7) is determined by the profiles of concentrations from 
the solution of BVP associated with the boundary condition (BVi), and Pi\x\ z\,d\) 
defined in (2.7) is determined by the profiles of concentrations from the solution of 
BVP associated with the boundary condition (BVo). 

Proof. For the experiment associated with the boundary conditions (BVi) where 
> 0 and = 0, formula (3.1) gives the influx of the tracer 


jM = jfltA e P [ i ] {p\z 1 ,d 1 ) 


i e p\ [x-,z\4i) 
0 £>i(x)/i(x) 


dx 


-i 


Similarly, for the experiment associated with the boundary conditions (BVo) 
where = 0 and R^’°^ > 0, formula (3.1) gives the efflux of the tracer 




1 ,di) 

0 T>i(x)/i(x) 


dx 


-l 


Formula (4.3) for the ratio of and Jf’ D| then follows 


rM 


□ 
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Corollary 4.2. Assume Hypotheses 2.1 and 2.2. For the one-isotope-setup (Setup 
1) with ideal ionic mixtures, the flux ratio of the two unidirectional fluxes is 

j[t,i] gZiCoVo rl e ziCofl°fx) , r 1 e ziCofl l fx) -l 

R lh°] J 0 Dffx)h(x) dx \J 0 Dfix)h{x) / 

where <(x) is t/ie profile of electric potential associated to the boundary condition 
(BVi) and (ft 0 \x) is that associated to the boundary condition (BVo). 

Note that, even for ideal ionic mixtures where (see (2.7)) 

p[\x-,zi,di) = and p^°\x\ z 1} d t ) = zffo^^x), 

the two profiles r/>W (x) and 4>^ (x) are different due to different boundary conditions 
in (BVi) and (BVo). This difference in the electrical potential reflects the effects of 
permanent charge, more than the difference in boundary conditions. Thus, even for 
ideal ionic mixtures, the PNP model shows that the flux ratio in the one-isotope-setup 
(Setup 1) is specific and contains information on the permanent charge. 

4.2 Flux ratio for the two-isotope-setup (Setup 2) 

For this setup, one uses one main ion species, one counter ion species, and TWO 
tracers T\ and T 2 ( two different isotopes of the main ion species ) in the SAME exper¬ 
iment. This is essentially the same setup used by Ussing in his derivation of formula 
(1.1). It turns out the flux ratio from this setup is universal when local models for 
excess potentials are used. 

The boundary concentrations for tracers are as follows: 

for tracer T\ : >0 at x = 0 and =0 at x = 1; 

for tracer T 2 : =0 at x = 0 and >0 at x = 1. 

This would yield an estimate for the influx of tracer T\ and an estimate for 
the efflux jf 2 ’°^ of tracer T 2 . We will examine the flux ratio 
In the following table, we summarize the notations used. 


Setup 2 

Concentration 

Valence 

Flux 

Main ion 

c i 

zi > 0 

tH 

j i 

Tracer T\ 

TPil 

C 1 

zi > 0 

ABA 

Tracer T 2 

r \h\ 

C 1 

zi > 0 

0*2,0] 

Counter ion 

C2 

Z2 < 0 

J-2 


The superscript [fi] in etc. indicates the corresponding quantity is 

associated with tracer T\, and the superscript i in indicates that it is the inward 
unidirectional flux of tracer T±. Similar remark applies to [ t 2 ] and o for tracer T 2 . 

Let (0(x), (x), ^(x), c^(x), c 2 (x), J^ n \ , J 2 ) be a solution of bound¬ 

ary value problem of system (2.1) with boundary conditions 

m = Vo, 4 ml (0) = L!, cf l] (0) = Lf l] , cf 2l (0) = 0, c 2 (0) = L 2 ; 

0(1) = 0, 4 rol (l) = Ri, cf l] (l) = 0, cf 2l (l ) = R l l 2 \ c 2 (l) = i? 2 . 
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The crucial observation is, since the two tracers are involved in the same experi¬ 
ment and they are assumed to have the same electrical and chemical properties (such 
as same valence z 1 and same diameter d\ , etc), both tracers have the same function 


pi{x\zi,di) = ziC,o(j){x) + 


1 

kflT 


H ex (x;zi,di) 


defined in (2.7) and the same diffusion coefficient D\. Concerning the flux ratio, we 
have the following. 


Theorem 4.3. Assume Hypotheses 2.1 and 2.2. For the two-isotope-setup (Setup 2) 
with two tracers T\ and T 2 of the same valence and same diameter in one experiment, 
the flux ratio between the influx jf 1 ’^ of tracer T\ and the efflux jf 2 '°^ of tracer T 2 is 


J 


[*!>*] 


L 


[h] o pi(0;zi4i) 


[*2 ,*] 


J 


gPi(l;zi,di) ' 


(4.5) 


Proof. Since the two tracers are involved in the same experiment and they have 
the same valence and diameter, the profile of the excess potential beyond the ideal 
component ksTpflx-, z\,d{) is the same for both tracers T\ and T 2 . Note also that 
the two tracers would have the same diffusion coefficient D\. 

Thus, for the tracer T\ with L ^ > 0 and = 0, formula (3.1) gives its influx 


j\P !>*] 


e pi(x;zi,di) 
o Di{x)h(x) 


dx 


-1 


Similarly, for the tracer T 2 with L ^ = 0 and > 0, its efflux is 


j|t2,o] _ _j^[t 2 \ e pi{i-,z 1 ,d 1 ) ^ J 


e pi(x;zi,dl) 


dx 


-1 


/ 0 Di(x)h(x) 

The formula (4.5) for the ratio of and then follows directly. 


□ 


Note that, with local models for excess components, pi(0; z\,d\) and pi(l; z\, d\) 
in formula (4.5) are determined by the boundary concentrations of the whole ionic 
solution, not just those of tracers. Nevertheless, when local models for excess compo¬ 
nents are used in PNP models, the flux ratio from this setup (Setup 2) is completely 
determined by the boundary conditions - independent of the channel structure. It is 
universal for all channel types (described by PNP type models) independent of the 
parameters of the channel. 


Corollary 4.4. Assume Hypotheses 2.1 and 2.2. For the two-isotope-setup (Setup 2) 
with two tracers in one experiment, if ionic mixture is treated as ideal ionic solution, 
then the flux ratio between the influx jf ] of tracer T\ and the efflux of tracer 

T 2 is 


J 


[*!.*] 


[* 2 , 0 ] 


r M 

,zT_gZiCoVo. 


J 


R 


[te] 


(4.6) 


Proof. For ideal ionic solutions, it follows from (2.7) that pi(x\zi,d\) = zi(of>(x). 
Formula (4.6) for ideal ionic mixtures is then a direct consequence of (4.5). □ 
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4.3 The one-isotope-setup (Setup 1) vs. the two-isotope-setup (Setup 
2) for flux ratios 

We comment more on differences among estimators of flux ratios from the two setups. 

At this moment, we are not able to give a comprehensive comparison between the 
flux ratios from the two setups. But, for the following situations, one may conclude 
that the flux ratio from the one-isotope-setup (Setup 1) is better than that from the 
two-isotope-setup (Setup 2) if the goal is to learn about the permanent charge. 

For simplicity, for the one-isotope-setup (Setup 1), we take, = p at x = 0 
and = 0 at x = 1 for (BVi) in (4.1), and l\' 0 ^ = 0 at x = 0 and R^'°^ = p at 
x = 1 for (BVo) in (4.2); for the two-isotope-setup (Setup 2), we take, = p and 
i)^ = 0 at x = 0, and = 0 and R = p at x = 1. 

(i) We already see that, for PNP with any local model for the excess potentials, the 
flux ratio from the two-isotope-setup (Setup 2) is given by (4.5) recast as follows: 

— pi(0;2i,<ii)-pi(l;zi,cZi) 

T [t 2 ,o] 

In particular, the two-isotope estimator has nothing to do with the channel structure. 
The two isotope estimator is universal. 

On the other hand, it follows from (4.3), the estimator of the flux ratio from the 
one-isotope-setup (Setup 1) is 


J 


[O] r , r 1 „p\ al (x;zi,di) 

1_ _ _ p v\ (0;zi,<b)-j4 ] (l;zi,di) / £_ 


•/ 


[ho] 

i 


= — e 


I o Di(x)h(x) 


dx 


1 e P\ (x;zi,di) . _i 

o Di{x)h(x) dX 


For PNP with local model for the excess potentials, we do have 
e p^(0;zi,di)-p( ol (l;zi,di) _ e pi(0;zi,di)-pi(l;zi,<ii) 


and hence, the first factor in the expression of from the one-isotope-setup 

(Setup 1) equals, itself, to the flux ratio f rQ m the two-isotope-setup 

(Setup 2). The last two factors in the expression of jM/jjbo] con ^ a j n significant 
information, namely, both factors involve permanent charge effects, even for PNP 
with ideal electrochemical components only. In this sense, the flux ratio measured in 
the one-isotope-setup (Setup 1) contains more information on channel structure than 
that in the two-isotope-setup (Setup 2). Although this came as a surprise to some of 
us, who had been inoculated with traditional views of unifluxes long ago ([21, 46]), it 
should not have. Setup 1 with one tracer involves two different spatial arrangements, 
and so it should give more information, particularly about any asymmetry in the 
charge distribution. 

Recall that p^(x;z\,di) is determined by the profiles of concentrations from the 
solution of BYP associated to the boundary condition (BVi) in (4.1) and (x; z\, d\) 
is determined by the profiles of concentrations from the solution of BVP associated 
to the boundary condition (BVo) in (4.2). The difference between and 

Pi\x;zi,d\) should be of order p. We are expecting that the channel structure 
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appears in terms of order p. This is indeed true as shown in formula (5.14) for the 
simple case that was treated in Section 5. 

(ii) As claimed in Remark 2.3, the above discussion is valid even if the local assumption 
(2.2) on the excess potentials is required only in the ionic solution outside the channel. 
We now consider PNP with nonlocal models for the excess potentials everywhere 
(assuming the BVP is well-defined). 

It is a reasonable assumption that the nonlocal model for the excess potentials and 
its local version differ in terms that are higher order (quadratic) in ionic diameters. 

Denote cq a characteristic concentration (number per length for one-dimensional 
model) of the ionic mixture. Let z be the valence and d be the diameter of the tracer. 

We will assume the ionic mixture is reasonably dilute near x = 0 and x = 1, and 
the permanent charge is small relative to the characteristic concentration. 

For clarity of the statement, denote the flux ratio from the one-isotope-setup 
(Setup 1) by ai(p,d,Q) and the flux ratio from the two-isotope-setup (Setup 2) by 
£72 (p, d, Q). Then, one can argue that, 

o l (p, d,Q) = — e z? ° y ° (l + ai(p, c 0 )(c 0 d) + a 2 (p, c 0 )(c 0 d) 2 + G(p,c 0 )— + h.o.t .), 
v cq J 

cr 2 (p, d,Q) = — e zCoV ° (l + ai (p, c 0 )(c 0 d) + a 2 (p, c 0 )(c 0 d) 2 + b(p, c 0 ) — ( c 0 d) 2 + h.o.t.). 

V cq J 

Note that the first three terms for both oq and cr 2 are the same and they DO NOT 
contain information of channel structure. Those terms are universal. The term 
G(p, co)^ for u i and the term b(p, co)^(cod) 2 for a 2 contain information of channel 
structure. They are specific. 

The dilute assumption would imply that (co<i) 2 <C 1 (for one-dimensional PNP). 
Therefore, permanent charge Q appears in the zeroth order 0(d°)-term in o\ but 
it appears in a 2 until the quadratic order 0(d 2 ). We conclude that the one-isotope- 
setup (Setup 1) is better than the two-isotope-setup (Setup 2) in capturing the effect 
of permanent charge in this situation. 

When the ionic mixture is crowded, the above formulas for a\ and o 2 are not valid 
and we would not be able to compare them in any definite terms at this moment. 


5 Case studies for flux ratios from the one-isotope-setup 


In this section, we consider specific flux ratios from the one-isotope-setup (Setup 1) 
in special cases to illustrate further some detailed information of channel structure 
that could be extracted from flux ratios. We assume 
(Al) the permanent charge is given by 


f 0 for x 6 (0, a) U (b, 1), 

1 Qo for x G (a, b ) 


where 0 < a < b < 1 and Qo is a constant. 

Note that the channel length has been normalized to 1 (see below for a dimen¬ 
sionless rescaling). 
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5.1 Rescaling of the model system (2.1) 

The main results for this special case rely heavily on the work in [51, 56] where 
the results are stated in terms of dimensionless variables after a rescaling. We will 
translate the results in terms of the original variables but for the reader’s convenience 
in comparing the statements, we present a full dimensionless rescaling below. 

The following rescaling (see [32]) or its variations have been widely used for con¬ 
venience of mathematical analysis. Set 

C 0 = max { sup |<2(x)|, L k , R k }, D 0 = sup {D k (x)}, i r = supe r (.x). 

x x,k x 


Let the channel length be l. Note that, in the previous sections, we have normal¬ 
ized the length of the channel to l = 1 and will still use this normalization except in 
the following rescaling. We make the re-scaling to get dimensionless variables 


2 i r £o k B T ~ D k (x ) 

6 = W ° k(x) = TC 

< i>(x) = = Co4>{x), 44) = Jk = 


A/ X Q{x) „ , . £r(x) 

Q(x) = £ r {x) = —— 

(_0 £r 

Jk 


(5.2) 


k B T TK ' ^ Co ' ” IC 0 D 0 

In terms of the new variables, the BVP (2.1) and (2.2) becomes 


d 


h(x) dx 

1 


i r (x)h(x 


d(j) 

dx 


= -^2z s c s - Q(x ), 
dJ ^ 


(5.3) 


k BT bk(x)hix)6td S = - Jk ' T = °’ 


with boundary conditions at x = 0 and x = 1 


m = coVb, 4 (o) = % m = o, 4(1) = 

C'O Lo 


(5.4) 


We will assume a special form for the permanent charge and include a hard-sphere 
potential. More precisely, we assume 

(A2) i r (x) and Dk(x) = D k are constants, and the dimensionless quantity £ is small. 
(A3) a one-dimensional local hard-sphere potential model is 

1 LHSf \ , C SC ] r ~ ( 4 _i_ dkY)j=i c oCj{ x ) 

hk (■ x) = - In \1 - 2^ djCoCjix)J + 

3 =1 


k B T 


1 - E?=i djCoCj(x) 


(5.5) 


+ dk)CoCj(x) + 0{d k dj ) 
j= 1 


where dj is the diameter of the jth ion species. 

The PNP model with local hard-sphere models (5.5) was studied in [54]. The 
local hard-sphere model (5.5) is an approximation of the nonlocal hard-sphere model 
adopted in [50, 58] from [73]. 
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The effects of the permanent charge and ion sizes on flux ratios will be examined. 
This will be accomplished by considering the permanent charge effect first and fol¬ 
lowed by ion size effect, and then the combination of these two. All these are done 
with extra assumptions so that a regular perturbation approach allows us to obtain 
detailed and useful information. Hopefully, this perturbation result provides insights 
and indications for more realistic models and for the real biological problem. While 
our analysis is neither general nor perfect, we (immodestly) believe it significantly 
more helpful than no discussion at all ([37, 77]) or analysis that assumes electrical 
potentials are independent of charge and structure. 

In the rest, we consider an ionic mixture of two (n = 2) ion species containing 
one main ion species with valence z\ > 0 and diameter d\ and a counter ion species 
with valence Z 2 < 0 and diameter da. Our assumptions are (A1)-(A3), in addition to 
Hypotheses 2.1 and 2.2. Set d = d\ and A = da/d\. For the local hard-sphere model 
(5.5) in (A3), in terms of the original variables, one has, from (2.7), 


pi{x; Zi, di) =z 1 Qo(t){x) + - 1 — Pi ns (x) 

k B T 

=zi(o<fi(x) + d(2c\{x) + (1 + A)c 2 (.t)) + 0(d 2 ). 


(5.6) 


5.2 Approximations of specific flux ratios 

Consider the one-isotope-setup (Setup 1) with one tracer but two sets of boundary 
conditions. For one experiment, the tracer is injected from the left boundary with 
concentrations = p and = 0 that produces the influx ^. In another 

experiment, the same tracer is injected from the right boundary with concentrations 
Rf’°^ = p and L^'°^ = 0 that produces the efflux . 

Our goal here is to provide an explicit approximation formula for the flux ratio. 
Of course, this is nearly impossible in general and only approximations are possible. 
Experimentalists are mostly concerned with particular cases, particular channels of 
biological and medical interest. Each of those cases will allow its own appropriate 
approximations because each case will be known in detail experimentally. Thus, the 
lack of generality is much less serious than it seems. The tools provided here will allow 
analysis appropriate for any particular system that is reasonably well studied. Our 
approximation has two levels. More precisely, with the assumption that e is small, 
we will only examine the zeroth order in e term of the flux ratio. (We will abuse the 
notion in the following to refer to the zeroth order in e term of the flux ratio simply 
as the flux ratio.) Then the flux ratio depends on ion diameters, permanent charge, 
and the injected concentrations of the tracer, and, with the assumption that these 
quantities are also small, we are able to obtain an approximate formula of the flux 
ratio with explicit leading order terms in these quantities. 

The following quantities will be needed to state the results. Recall that h(x) is the 
cross-section area of the channel over x, and Q(x) = Qq for x 6 (a, b ) and Q(x) = 0 
for x 0 [a,b\. We denote H(x) = JJ h^ 1 (s)ds. For ease of notation, we set 


s 


Li _ H{a) 

/>v a ~W)' 


p 


H(b) 

H(iy 


(5.7) 
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and define fi(s), f 2 (V 0 ,s), 31 (Vo, s, a, /3), and g 2 (s,a,/3 ) as 


/i(») 


2 s In s — s 2 + 1 
s(s — 1 ) In s ’ 


f 2 (V 0 ,s) 


s + 1 

s(^iCoVb + Ins) 


gZiCoVo _|_ j 
se^iCo Vo _ i ’ 


(5.8) 


gi(V 0 ,s,a,/3) = - 


(/? — a) [(a + P — 2 a/3)(s — l ) 2 + 2 s] 


(21 - 22)[(1 - a)s + a] 2 [(l - p)s + /3] 2 (lns ) 2 

_ 2 (/?-a)(s 2 -l) _ 

{zi - z 2 )[(1 - a)s + a][(l - P)s + /3]s(lns ) 3 
(/3-a)(s + 1 ) 


{z 2 Co Vo 


(%2 Co Vo + 
Z2 Co Vo 


(21 - 22 )[(1 - a)s + a][(l - p)s + /3](s - 1) Ins 
(2 In s + s — 1 /s) [ ln((l — p)s + p) — ln((l — a)s + a)] 


(21 - 2 2 )(s - l) 2 (lns ) 2 

( s + 1 (1 — P)s + P 

32 (s, a, p) =7 - —7 - -r -—In -—— 

(zi — Z 2 )s(s — 1 ) In s (1 — ajs + a 

{P ~ a)(s 2 — 1 ) 


- In ,■ 

In s 
2 


22 Co Vo, 


(21 - 2 2 )[(1 - a)s + a][(l - p)s + /3]s(lns ) 2 ' 


(5.9) 


Denote the boundary concentrations of the main ion species and the counter ion 
species, before the tracer is added , by (. L\,L 2 ) at x = 0 and (Ri,R 2 ) at x = 1, and 
assume the electroneutrality conditions among them 

z\L\ + z 2 L 2 = 0 and z\R\ + z 2 R 2 = 0. (5.10) 

For the electroneutrality boundary conditions, we consider the following two cases. 


5.2.1 Total electroneutrality boundary conditions 


The first case is, after the tracer is included , one requires the total electroneutrality 
boundary conditions among all three ion species; that is, the boundary conditions for 
the two experiments are, respectively, 


(BVi) : 


' 0 [i] (O) = V 0 , 

c [ r\ 0)= Li t 

cf’ i] (0) =p>0, 

. <f(0) = l 2 - %p, 


</» [il (l) = 0 , 

c [ r\±)=Ri, 

cf’ i] (l) = 0, 

cf(i ) = R 2 ; 


(main ions) 

(tracer) 

(counterions) 


(5.11) 


(BVo) : 


' 0 [O] (O) = V 0 , 

c [ r\ 0 ) = L\, 

cf’ o] (0) = 0, 

. 4 o 1 ( 0 ) = L 2 , 


0 [o] (l) = O, 

c [ r o] a) = Ri, 

cf’ o] (l)=p>0 , 

C [° ] (l ) = R 2 -f 2 p. 


(main ions) 

(tracer) 

(counterions) 


(5.12) 


That is, for (BVi), one adds the extra amount —zip/z 2 of counter ion concentration 
at x = 0 to compensate the inclusion of the tracer concentration p at x = 0 ; and, for 
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(BVo), one adds the extra amount —zip/z 2 of counter ion concentration at x = 1 to 
compensate the inclusion of the tracer concentration p at x = 1. Therefore, 

For (BVi): zflLi + p) + z 2 (L 2 - zip/zf) = 0 and z\R\ + z 2 R 2 = 0, 

For (BVo): z\L\ + z 2 L 2 = 0 and zflRi + p) + z 2 {R 2 - z\p/z 2 ) = 0. 

Then we have 


Theorem 5.1. Assume d, p, and |Qo| are small. The flux ratio between the influx 
j\ L ’^ associated to the boundary condition (5.11) and the efflux associated to the 
boundary condition (5.12) is, with s = Li/Ri, 

tM 
■ J \ _ 

tM 

- e^° y ° (h(s) + / 2 (Vb,s)) - e^gflVo, s, 

Til it] XL 1 

+ 0(p 2 ,pd, d 2 , pQl , p 2 Qo,pdQ 0 ), 

where a and (3 are in (5.7), f\ and f 2 are in (5.8), and g\ are in (5.9). 

Formula (5.14) will be established in the Appendix (Section 7). The flux ratio 
is specific due to the third term involving the key parameters Qo, a and fl of the 
prescribed channel structure (see (5.1) and (5.7)). 


e ziCoVo ^ + ( 2(Ll _ Rl) + (i + A)(L 2 - R 2 ))d) 


5.2.2 Partial electroneutrality boundary conditions 


The other case is to simply keep the concentrations of the main ion species and the 
counter ion species for both experiments after the tracer is included; that is, the 
boundary conditions are 


(BVi) : 


</> [i] (0) = Vo, 

4^(0) = ^, 

Jt,i 


[t,i] 


cr(o)= P >o, c r j (i)=o, 

c5 ;1 (0) = L 2 > (1, c\\l) = R 2 >t). 


(mam lonsj 

(tracer) 

(counterions) 


(BVo) 


<^(0) = Vo, 


4 m,o] (0) = Li > 0, c‘ m ’ 0| (l) = Ri > 0, 


] ( 1 ) = 0 , 


m,o] , 


c i’° l (o) = 0, 


[ c'° ] (0) = L 2 > 0, 
In this case, one has 


cf’ OJ (l) = p>0, 

c'° ] (l) = R 2 > 0. 


(main ions) 

(tracer) 

(counterions) 


For (BVi): z\(L\ + p) + z 2 L 2 = zip « 0 and 2 :ii?i + z 2 R 2 = 0, 
For (BVo): ^iLi + z 2 L- 2 = 0 and z\(R\ + p) + z 2 R 2 = zip ~ 0. 


(5.15) 


(5.16) 


(5.17) 


Thus, one only has approximate electroneutrality boundary conditions after the tracer 
is added. But, since the amounts of tracers injected are small and the net amounts 
of charge resulting from approximate electroneutrality are exceedingly small, this 
approximate electroneutrality boundary condition is realistic and more practical too. 
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Theorem 5.2. Assume d, p, and |Qo| are small. The flux ratio between the influx 
associated to the boundary condition (5.15) and the efflux jf'°^ associated to the 
boundary condition (5.16) is, with s = L\/R\, 


M 

l_ 

M 


_ e ziCoVo ^ + ( 2 ( Ll _ Rl) + (! + a)(L 2 - i2a))d) 

_ e ziCoVo (^^ fl{s) + h{ y^ s) \ _E_ 

\Zi — Z2 J Hi 

- e zlCoKl — (gi(V 0 ,s,a,P) + 3 2 (s, 

+ 0(p 2 , pd, d 2 , pQl,p 2 Qo , pdQo), 


(5.18) 


where a and fl are in (5.7), fi and fi are in (5.8), and gi and g 2 are in (5.9). 

The derivation of formula (5.18) will also be given in the Appendix (Section 7). 


5.3 Comments on formulas (5.14) and (5.18) 

The derivation of formulas (5.14) and (5.18) in Section 7 is rather complicated. To 
gain additional confidence in the correctness of the formula, we have cross-checked 
it using a symmetry of the problem, that is, the result should be symmetric if one 
swaps the channel between the left end and the right end. More precisely, if we set 
x = 1 — x, then x = 0 becomes x = 1, and x = 1 becomes x = 0. If we denote the 
relevant quantities for the new situation with an overhat, then Lj = Rj, Rj = Lj, 
V 0 = —Vo, s = Li/Ri = 1/s, and etc.. Therefore, one 

should have 


fM tM ' 

■ J \ 

The latter is equivalent to, up to the leading orders in (5.14) and (5.18), the relations 


/i(l/s) = -sfi(s), f 2 (-V 0 ,l/s) = s /2 (Vo, s), 

gi(-V 0 ,l/s,l - fl, 1 - a) = -s 2 gi{Vo,s,a,fl), 

32 (l/s, 1 - fl, 1 - a) = - s 2 g 2 (s , a, fl) 

where f\, f 2 , gi and 52 are defined in (5.8) and (5.9). We have verified them all! 

For ideal ionic solutions, formulas (5.14) and (5.18) reduce to, for the total elec¬ 
troneutrality boundary conditions (5.11) and (5.12), 

(/,(*) + A(V4,.) + »(V4,.,«,/3)|) £ (519) 

+ O(p 2 ,pQ 2 0 ), 


and, for the approximate electroneutrality boundary conditions (5.15) and (5.16), 


J 


M 


1 — _ gZiCoVo _ gZiCoVo 


J 


[ho] 


^/l(5) + / 2 (Fo,s)j£ 


— e 


Cob,-— (gflVo,s,a,fl) + g 2 (s,a, + 0(p 2 ,pQl) 

Z\ — Z2 \ /ill Ra 


(5.20) 


R\ R.\ 
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It is noted that, even for ideal ionic solutions, the flux ratio from the one-isotope- 
setup (Setup 1) is specific and captures some properties of the channel structure. 

Implications of formula (5.14). The following properties of /i, fi, g\, and 52 can 
be established easily and will be used in the discussions below. 


Lemma 5.3. For 0 < s < 1 , fi(s) > 0, for s > 1 , fi{s) < 0, 


lim fi(s) = + 00 , 

s— >0+ 

lim <71 (Vo, s, a, f$) 

s—>1 


lim 92 (s, a, /3) 

s —>1 


lim /i(s) = lim /i(s) = 0, _ lim f 2 (V 0 ,s ) = 

s —>-1 s—too 21 Co Vo—^—Ins ZS 

- ((2a + 2/3 ^ 1)2 + ^ ^ - !) 2 - !]) 

2(/3 — a) (a + j3 — 1 ) 

Zl - ^2 

_ (/3-a)(q + /3- 1) 


We discuss some indications of each of the three leading terms on the right-hand 
side of the formula (5.14). Similar comments apply to formula (5.18). 

(a) The first term 


_ e *iCob> (1 + ( 2(Ll _ Rl) + (i + A)(L 2 - i? 2 ))d) 


represents directly the asymmetry of the electrochemical potentials (the concentration 
and the electrical ‘driving force’ using the language of the laboratory, introduced by 
Hodgkin and Huxley) at the two boundary points x = 0 and x = 1 (independent of 
the concentration p of the tracer). In particular, as p —> 0 in (5.14), 


tM 

lim - e Zl(oVo - e Zl ^° v ° (2(L X - R^) + (1 + A)(L 2 - R 2 ))d + 0{d 2 ). 

p^r 0 Jp°\ 

Note that, the value for the limit of the flux ratio, NOT equaling to —1 in general, 
reflects the asymmetry of the boundary conditions. The effects caused by permanent 
charge disappears in the p —> 0 limit. A direct justification (as an alternative to a 
consequence of the formula (5.14)) of the latter is given in Remark 7.1. Remark 7.1 
also explains why there is no dQ 0 term in the approximation formula (5.14). 

(b) In formula (5.14), the term 


-e*i^o(/i( s ) + / 2(K 0 , s ))-|- 

represents the interaction between the tracer and the asymmetry of the boundary 
electrochemical potentials. 

(c) In formula (5.14), the term 

— e Zl ^ v ° gi (Vo,s,a,P)^ 

W 

contains the interaction between the tracer, the asymmetry of boundary electrochem¬ 
ical potentials, and the permanent charge (in this simplified setting). 
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(d) If the boundary conditions are symmetric , that is, Vo = 0, L\ = R\ and L2 = R2 
then, from Lemma 5.3, formula (5.14) reduces to 


J 


[t,i\ 


J 


[ho] 


= -1 + 


2(/3 - a)(a + /3 - 1) Q 0 p 
z 1 - Z2 R\ R\ 


+ 0(p 2 , pQo)- 


In this case, the effect of permanent charge on the flux ratio is particularly captured. 

Flux ratio and the flux ratio exponent n'. Next, we discuss the implications of 
formula (5.14) for a relation between the flux ratio with the flux ratio exponent n! 
introduced by Hodgkin and Kenyes in (1.3), under further assumptions. Similar 
comment applies to the term in formula (5.18). 

Recall, from (1.3) with = Rf['°^ = p, that 


■[Pi] 

1_ _ _ p ziC 0 V 0 n' 

M 


Therefore, the flux ratio exponent n! is given by, ignoring the higher order terms, 


1 


n 


zi( 0 Vo 


In — 


J 


[t,i] 


J 


[ho] 


T + 
+ 


g 1 (Vo,s,a,/ 3 )Q 0 /Ri + fi(s) + f2{V 0 ,s) p 


Z 1 C 0 V 0 

2 d\(L\ — R\) + (di + d2){L2 — R2) 
^iCoVq 


Ri 


In particular, if the concentrations Lf s and Rf s at the boundaries are such that 
2di(Li — R\) + {d\ + d2){L2 — R2) can be ignored, then one has 


/ - gi(Vo,s,a,/ 3 )Qo/Ri + /i(s) + jgQjh 5 ) P 

ziCoVo R\ 

Since the above formula is obtained under the assumption that Q 0 is small and 
the ionic mixture is reasonably dilute, we cannot apply it to large Qq. BUT, if we 
pretend that the formula were true for large Qo, then 51 (Vo, s, a, fi)Qo/R\ dominates 
fi(s) + / 2 (Vq, s), and hence, 


, , . gi(V 0 ,s,a,/ 3 )Q 0 p 

n ~ 1 + —pm, —a?' 


which implies, when g i(Vq, s, a, (3) / 0, 


/ . , r gi{Vo,s,a,P) 

n > 1 if -—. -Qo > 0 

ZiCoVo 


and n < 1 


gi (y 0 ,s,a,P) 

'if --Qo < 0; 

ziCoVo 


that is, n! > 1 for Qo with one sign and n! < 1 for Qo with the opposite sign. 

Detecting (a,( 3 ,Q 0 ) using ( 5 . 14 ). Assume the permanent charge has the general 
structure characterized by (a,/3,Q 0 ). We will propose a procedure for determining 
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these quantities. Rearrange formula (5.14) to get 

+ l+(2(Lx - R 1 ) + (1 + A)(L 2 - R 2 ))d + (/i(s) + f 2 (V 0 , s))^- 

J x ’ i 

= - c/i(R 0 ,s,a,/3)^y2 +0(p 2 ,ph, d 2 , pQ 2 0 , p 2 Q 0l pdQ 0 ). 

All terms on the left-hand side can be determined either by experimental setup or 
by measurement. Under the assumption, the first term on the right-hand side is the 
leading term that contains information on net or overall effect of permanent charge 
and channel geometry through its dependence on Q o and (a, j3). In general, three 
choices of (Vo, L\, R\) lead to three equations for (Qo, a, f3) that could determine the 
values of Qo, & and (3. 

6 Discussion 

As an extremely important experimental method in studying of channel properties, 
flux ratios have been widely applied ([37, 77]). Their dependence on physical param¬ 
eters is subtle and is dictated by details of experimental designs. To emphasize the 
importance of the subtleness of flux ratios, we classify them into two main classes: 
universal and specific. Universal flux ratios are those independent of channel struc¬ 
tures and specific flux ratios contain information of channel structures. 

In this paper, we have examined two applications of flux ratios: (i) proportionality 
constant of the flux of a main ion species and the flux of its tracer; (ii) flux ratios of 
tracers associated to two experimental setups. 

For the first topic, it is shown that, based on analysis of PNP type models in¬ 
cluding local excess potentials in Section 3, the proportionality constant of the flux 
of a main ion species and the flux of its tracer is universal and can be completely 
determined by boundary conditions, and hence, one can estimate the flux of the main 
ion species from that of its tracer. Furthermore, an experimental design is proposed 
(Remark 3.3) based on the analysis so that the proportionality constant can be de¬ 
termined in a simple and precise way even for PNP models with nonlocal models of 
excess potentials. 

For the second topic, based on analysis of PNP type models with local excess 
potentials in Section 4, we examined two popular experimental designs for measuring 
flux ratios of tracers. We showed that the two-isotope-setup (Setup 2) yields universal 
flux ratios that do not reveal channel structure but the one-isotope-setup (Setup 1) is 
specific even with dilute ionic mixtures. As claimed in Remark 2.3, it is not hard to 
see that the local nature in Hypothesis 2.2 can be replaced by assuming the hypothesis 
only at the bathes as long as the BYP is well-posed. 

Presumably, if one takes PNP with nonlocal excess potential models (inside and 
outside of the channel), then the flux ratios for both setups might capture information 
on channel structures. However, at this moment, it is not clear how different they 
are, and this requires a serious investigation of PNP systems with nonlocal excess 
potential models. To our best knowledge, PNP including nonlocal excess potential 
models has not been mathematically analyzed. In fact, the well-posedness issue for 
such a model is not clear due to its infinite number of freedoms; that is, what should 
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be “the boundary conditions” in order for the boundary value problem of a PNP 
model with nonlocal excess potentials to have a unique or finite many solutions? An 
initial effort is taken in [75] for a simplified setup and hope this activity can stimulate 
further investigations on PNP model with nonlocal excess potentials. 


7 Appendix: Derivations of results in Section 5.2 

In this section, we provide a derivation of formula (5.14) in Theorem 5.1 and formula 
(5.18) in Theorem 5.2 from the basic formula (4.3) in Theorem 4.1 for flux ratios from 
the one-isotope-setup (Setup 1). The derivation for both formulas are similar so we 
will provide the detailed derivation for formula (5.14) in Section 7.2 and comment on 
differences in the derivation for formula (5.18) in Section 7.3. 


7.1 A reformulation of (4.3) for cases in Section 5.2 


We will first recall the setup of the case specified in Section 5.2 and give a new form 
of formula (4.3) for the special case, which is convenient for us to apply previously 
established results in late part. 

Recall that, for the one-isotope-setup (Setup 1), one tracer is used for two exper¬ 
imental settings with two different sets of boundary conditions for tracers. For one 

\t i] 

experiment, the tracer is injected from the left boundary with concentration L l f = p 
and zero right boundary concentration = 0 that produces the influx In 

another experiment, the same tracer is injected from the right boundary with con¬ 
centration = p and zero left boundary concentration L ^ = 0 that produces 

the efflux jf'°\ 

Formula (4.3) in Theorem 4.1 for the flux ratio between jj f '^ and jf' 0 ' 1 is 


J 


M £tM] e j4’l(0;zi,di) rl e p( o] (x;zi,di) 


J 


[ho] 


dx 


R [t’°] e p l ° ] (i- )Zl ,di) Jo Di(x)h(x) \Jo Di(x)h(x) 


"! pPi'{x\zi,<h) \ —l 

-dx] , (7.1) 


where p 1 ^ (x; z\,d\) defined in (2.7) is determined by the profiles of concentrations from 
the solution of BVP associated with the boundary condition (BVi), and p^°\x; z\, d\) 
defined in (2.7) is determined by the profiles of concentrations from the solution of 
BVP associated with the boundary condition (BVo). 

Remark 7.1. Note that, as p —> 0, one has (x; z\,d\) = p^°\x; z\, d\), and hence, 


lim 


J 


[t,i\ 


P^o jV 


[ho] 


(0;zi,di) 


in particular, this limit does not contain information on permanent charge Q. An 
technical implication is that all terms involving Qo in the approximation of 
should have a factor of p. This explains why there is no dQo-term in the leading order 
expansions 0 /j|h*l / j|*>°] formula (5.14) in Theorem 5.1. 
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For the case studied in Section 5.2, we are able to obtain an approximation for the 
flux ratio, up to leading orders, explicitly in terms of given quantities of the problem. 
It follows from (3.1) in Proposition 3.1 that 


O fl e p [ i ] (x\zi,<h) , _i 

o Di(x)h(x)J 


J' 


(L\ + Lj’^)e p i J(0;zi,dl) 


i2 1 e p i 1 (i; z i ,<i i) 5 


where ./[^ = + jf' 1 ^ is the sum of the fluxes of the main ion species and its 

tracer. Similarly, 


r-i gP^fxjZi.di) Lie pl ° ] ( 

-dx = ---n-----. 


/o Di{x)h{x) 


J 


l 


where is the sum of the fluxes of the main ion species and its 

tracer. Therefore, 


tM 

yi_ 

7-[ho] 

J 1 


-(/) • (//) • 


«/' 




(7.2) 


where the factors (I) and (II) are given by 


(I) 


^M e p( ll (0;zi,di) 


(II) 


Lie Pi ] (°-’Zi,di) _ (i?! + i?M)ePi oI ( 1 ^i> d i) 

(Li + Lf’*]) e Pi ] ( 0 ^i> d i) - iijePi 1 ! 1 ;*!’*) 


(7.3) 


The factors (I) and (//) on the right-hand side of (7.2) are determined by bound¬ 
ary conditions and are independent of the permanent charge. We will first determine 
these two factors and determine the last factor ./[^/,/[ 0 ^ afterwards. 


7.2 Derivation of formula (5.14) 

For this formula, the electroneutrality boundary conditions (5.13) are assumed to 
hold among all ion species including the tracer. It follows from (5.6) that, for the 
boundary condition associated to (BVi) in (5.13), 

p^(0;zi,di) =zi( 0 V 0 + d{2(Li + p) + (1 + A )(L 2 - z 1 p/z 2 )J + 0(d 2 ), 
=~iCofo + d(2L\ + (1 + A)L 2 ) + 0(d 2 , dp), 
pf(l;zi,di) =d(2R 1 + (1 + \)R 2 ^j + 0(d 2 )- 

and, for the boundary condition associated to (BVo) in (5.13), 

p\ ^ (0; zi, d\) =^iCol / o + d(2L\ + (1 + \)L- 2 ) + 0(d 2 ), 
p [ °\l;zi,di) =d(2(Ri + p) + (1 + A )(R 2 - zip/z 2 )) + 0(d 2 ) 

=d(2Ri + (1 + X)R 2 ) + 0(d 2 , dp). 
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Thus, the factor (/) on the right-hand side of (7.2) is 


(/) =e p i 1 ( 0;Z1 > dl )-Pi’ 1 (Mi > d i) 

=e*i& Vo ( 1 + (2(L 1 -i2 1 ) + (l + A)(L 2 -i2 2 ))d) + 0(d 2 ,dp), 
and the factor (II) on the right-hand side of (7.2) is, with s = L\/R \, 


(7.4) 


L\e Zl ^° v ° - (i?i + p) + [(2L 2 + (1 + A )L 1 L 2 )e z ^° v ° - (2R 2 + (1 + A)i?ii7 2 )] d 

(Lx + p)e z iCoVb - + [(2Lf + (1 + A )L 1 L 2 )e z R°v 0 _ ( 2 j r2 + (i + A^ifc)] d 

+ 0(p 2 ,pd,d 2 ) 


e ziCoVo + 1 p 
se ziC, 0 V 0 _ i f> A 


+ 0(p 2 ,pd,d 2 ). 


(7.5) 


It remains to approximate the factor ./[^/in (7.2) that depends on the per¬ 
manent charge Q and hence the full profile of the electric potential and the ion con¬ 
centrations. Its approximation is complicated even for the simple case we considered 
in Section 5. We start with a general expression 


A 1 .1 

= S Q (p) + Si(p)d + S 2 (p)Qo + 0(d 2 , pdQo, pQl), 

and determine the quantities Sk(p )'s for k = 0,1, 2. In particular, we are interested 
in So(p ) up to O(p) order, S\(p) = Si(0) + O(p) due to the factor d in S\(p)d, and 
S 2 (p) up to O(p). The reason for the latter is, from Remark 7.1, any term involving 
Qo should have a factor p. The advantage of this expression is as follows. One can 
assume the ionic solution is ideal to get the term So(p ) + S 2 (p)Qo (Section 7.2.1) 
and then assume the ionic solution is nonideal but the permanent charge is zero to 
get S 0 (p) + S\(p)d (Section 7.2.2). The superposition of them then provides the 
approximation for ./[^ /j\°^. 


7.2.1 Case with ideal ionic solution for So(p ) + S 2 (p)Qo 

We will apply results from [51] for the classical PNP for this case. 

It follows from results in [51] (displays (4.4) and (4.5) in [51]) that the influx of 
the tracer for the boundary concentrations in (BVi) is 

j¥ ] =J® + J$Q 0 + O(Q 2 ), (7.6) 

where 

j[i ] _ D\(L\ + p — R\) Mi’* 

10 ~H(l)(ln(L 1 + p)-lnR 1 )k B T' 

j[i\ _ DiA^(z 2 (l - B®)CqVq + ln(Lj + p) - In Rj) p, 6 / (7 ' 7) 

(zi - z 2 )H(l)(h\(Li + p) - InRi) 2 ^bT 

and 

i‘ “ A^k 1 )) = ^iCoVq + ln(Ti + p) - InRi, (7.8) 
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is the difference between the boundary electrochemical potentials for (BVi), and 

A [i] = _ (P ~ a )(Li + P — Ri) 2 _ 

((1 — a)(L\ + p) + ai?i)((l — P)(Li + p) + pRi) In ~, 

p>[i] = _}_ . (1 - P)(Li + p) + {3Ri 
(1 — a)(L\ + p) + aR\ 


Similarly, the efflux of the tracer for the boundary concentrations in (BVo) is 


JP = 


=jS + j{i ] Qo + C>(^), 


(7.10) 


where 

j[o] Di(Li~(Ri+p)) 

10 //(1)(In /.! — hl(//] -p))k,ir 

j[o] _ DiAW(z 2 (1 - BM)C qVq + lnh, - ln(gi + p)) ^° (7 ' 11} 

11 (z 1 -z 2 )H(l)(lnL 1 -ln(R 1 + p)) 2 k B T ’ 

and 

^/4’° : = ^ (> ] (0) - h [ ° ] (l)) = ziCoVo + hr L x - ln(Ri + p) (7.12) 
is the differences between the boundary electrochemical potentials (BVo), and 


A [o] = _ (P-a)(Li - Ri ~ p) 2 _ 

((1 — a)L\ + a(R\ + p))((l — P)L\ + P{R\ + p)) In 

D[o] _ 1 , (1 - P)L\ + P(Ri + p) 

A& (1 — a)L\ + a(i?i + p) 

Therefore, 


jf jg + 2 _ ( , n 

T “7H + jHq 0 + 0(<3o) 0 V + + 0(Qo) - 

where 

_ _ Li + p - Ri In Li - ln(Ri + p) /x^ 

° J io Ll ~ ( Rl + P ^ ln ^ Ll + P)~ ln Rl l4'° ’ 

2^ jW jH ^ 2(1 — 7?W)CoVo + ln(Ti + p) — lni?i^ 

F ° J io J Io ( Zl ~ Z2 ^ Ll + P~ Ri H ln ( Ll + p)~ lni?i) 

A® (z 2 ( 1 - B[°])CoV 0 + In- ln(i?! + p)) 

(zi - z 2 ){Li - (Ri + p))(lnLi - ln(i?i + p)) 
Lemma 7.2. With s = L\jR\, one has 


F 0 =l+[fi(s) + 


8 + 1 


s(ziCoVo + Ins)/ Ri 


+ 0(p 2 


(7.13) 


(7.14) 


(7.15) 
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where fi(s) is defined in (5.8), and 


It =9i (Vo,s,a,f3)-^ + 0(p 2 ), 
ro iiI 

where gi(Vo, s, a, j3) is given in (5.9). As a consequence, 

S«(p) = l+(h (») + s( 2 lCo V^+ lns) )£ and SM= gi (V 0 , S , a ,n^. (7.16) 

Proof. For the term Fq, it follows (7.15), (7.8) and (7.12) that 

_s-l + lns-ln(l + ^-) ziCoVo + ln(s + ^-) 

° s - 1 - ^ In s + ln(l + j|-) Co Vo + In s - ln(l + ’ 

_ /2slns - s 2 + 1 s + 1 \ p . 2 \ 

\ s(s — 1) In s s(ziCoVb + Ins)/ R\ P 

which is the claimed formula for Fq. 

The approximate formula for F\/Fq is lengthy but otherwise straightforward. 
With s = L\/R\ and = Rf’ 0 ^ = p, it follows from (7.15) that 

F\ _ z 2 (AW-AWbM )CoVo AM 

Fq {z\ - z 2 )R\(s - 1 + p/R\)iy In s + p/(sR \)) + (z\ - z 2 )R 1 (s - 1 + p/R\) 

_ Z2(AW - AWbW)C qVq _A|°l_ 

{z\ - z 2 )R\(s - 1 - p/i?i)(lns - p/Ri) (z\ - z 2 )Ri{s - 1 - p/R\) P 


Next, we fix A W, B^\ Al°l and B^°\ and expand the above in p to get 


F\ 

Fq 


z 2 (qVq{AW - - AM + AMsW) (s - 1)(AW - AH) 

(z\ - z 2 )Ri(.s - 1) In s + (zi - z 2 )Ri(s - l) 2 

Z 2 <ZoV 0 [(AW - AWsW)(lns + 1 - 1/s) + (AH - A^B^)(lns + s - 1)] p 

(zi - z 2 )(s - l) 2 (lns) 2 R\ 


AW + AH p 
(zi - z 2 )(s - l) 2 R\ 


+ 0(p 2 ). 


Now, we expand AW, A^°\ £>W an d #[°] dehned in (7.9) and (7.13) in p as 


AW =aJ 1 + Af ] + 0(p 2 ), AM = A [ ° ] + A [ f + 0(p 2 ), 

rtl -tii 

A®B® =(A®B% + (AWsW)i-|- + 0(p 2 ), 

Tti 

aH^H = (aH j b[° 1 )o + (AWs[°])i-|- + 0(p 2 ). 
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It follows directly that 


4 1 =4° ] = -(P-a) 

=-(/?- a) 

+ (/? — a) 

+ (/? — a) 

4°’ =(P ~ a) 


(s-1) 2 


[(1 — a)s + a} [(1 — f3)s + /3\ In s ’ 

_ 2(«-l) _ 

[(1 — a)s + a][(l — /3)s + (5} Ins 
(s — 1) 2 [2(1 — a)(l — f3)s + a + /? — 2 a/3\ 
[(1 — a)s + a] 2 [(l — f$)s + (3} 2 Ins 

_ (^-l) 2 _ 

[(1 — a)s + a][(l — /3)s + /3]s(lns) 2 ’ 
2(s-l) 


+ ((5 — a) 

~ (P ~ «) 


[(1 — a)s + a] [(1 — (3)s + (3\ In s 
(s — l) 2 [(a + f3 — 2 a/5)s + 2 a/3] 


[(1 — a)s + a] 2 [(l — /3)s + /3] 2 In s 
(s-1) 2 


[(1 - a)s + a][(l - f3)t + /3](lns) 2 : 

(aWb'-])o =(A[°laM) 0 ^ In 

(1 — Glj<S + CX 

(Al‘lBW)i l ~ f} 1_ “ 


(1 — P)s + (3 (1 — a)s + a’ 


(A^bM) 1 


p 


a 


(1 —/3)s + f3 (1 — a)s + a 

In particular, 

Af - A [ ° ] =0, (A®BV ] )o - (A^B^)q = 0, 


(7.17) 


4? - (A®B®) o = - 

4> - TlP = - 
+ 

+ 


(/3 - a)(s - if 


[(1 — a)s + a][(l — /3)s + f3\ Ins 
4(/3 — a)(s — 1) 


— In 


(1 ~ P)s + P 

(1 — a)s + a 


[(1 — a)s + or][(1 — /3)s + 0\ Ins 
(/3 — a)(s — 1) 2 [(2 — 3a — 3/3 + 4 a/3)s + a + (3 — 4 a/3] 
[(1 — a)s + a] 2 [(l — /5)s + /3] 2 Ins 
(P - a)(s + l)(s - l) 2 


[(1 — a)s + a][(l — j3)s + /3]-s(lns) 2 ’ 

<- 4l,lBl "> 1 - < aIoIb,o1) ‘ = 0^4 - 44 W 
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Combining the above estimates, one has 

F\ _ _ 4(/? - a)z2CoVo _ P_ 

F 0 {zi - z 2 )[( 1 - a)s + a][(l - f3)s + /3](lns) 2 (Ri) 2 

(f3 — a)z 2 (oVo(s — 1)[(2 — 3a — 3/1 + 4 a(3)s + a + f3 — 4 a/3] p 
(zi - 2 2 )[(1 - a)s + a] 2 [(l - f3)s + /?] 2 (lns) 2 (Ri) 2 

_ (ff - a)z 2 (oVo(s + l)(s - 1) _ p_ 

(zi - z 2 )[(1 - a)s + a][(l - /3)s + /?]s(lns) 3 (Ri) 2 

_ ^ r 1-2/3 _ 1 — 2a i p 

(zi — z 2 )(s — 1) Ins .(1 — j3)s + /3 (1 — a)s + a. (Ri) 2 

_ (fl ~ a)z2(oVo(21ns + s - 1/s) _ 

{z\ - z 2 )[(l - a)s + a][(l - (3)s + /3](lns) 3 (Ri) 2 
Z2CqVq{2 In s + s - 1/s) (1 -/3)s + f3 p 

(z\ — z 2 )(s — l) 2 (lns) 2 (1 — a)s + a (Ri) 2 

(/3 — a)[(a + (3 — 2 a/3)(s — l) 2 + 2s] p 
{zi - z 2 )[( 1 - a)s + a] 2 [(l - (3)s + /?] 2 lns (Ri) 2 

,_ — a)(s + l)(s — 1) _ p 2 

(z\ - z 2 )[{ 1 - a)s + a][(l - (3)s + /?]s(lns) 2 (Ri) 2 

The fifth term above can be split and combine with the first and the third terms 
to get the formula for F\/F$ claimed in the statement of the lemma. □ 


7.2.2 A case with a hard-sphere component for So(p) + S\(p)d 

In this section, we examine the flux ratio for PNP models with a hard-sphere com¬ 
ponent (point-charge with volume exclusion) and assume zero permanent charge. 
Recall the local hard-sphere potential (5.5) 

1 LHS / \ , , / A , dkYJj=i c j( x ) 


(X) = " 111 l 1 - E d &W) + , >:; { <l i c i (.r) 


= ^{dj + dk)cj(x) + 0(d k dj ) 

3 =1 

where dj is the diameter of the jth ion species. 

For two ion species with diameters d\ and d 2 , set d = d\ and d 2 = A d, and expand 
the fluxes as J k = J k Q + J k \d+0{d 2 ), k = 1, 2. As observed in [29], the Nernst-Planck 
equations imply that the flux J* is proportional to the transmembrane electrochemical 
potential pf = Pi(0) — Pi( 1) where 

-'—pi =ziCoV& + hr Lx - In R\ + ~ (1 + A)d ^ + 

ktsT z 2 

7—tt;/4 =z 2(oVq + lnL 2 — lnR 2 — --—- ~{L 2 ~ R 2 )d + 0(d 2 ). 


(7.18) 


It follows from results in [50, 54], for example, from Corollary 3.6 in [54] (with 
different but equivalent expressions) that the flux ./[^ is given by 

j? = (1 + 2 F^1 W f'h + o(cP), 

V z\z 2 J ziH(l) k B T 
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where, with p = = Ri’°\ 

-—— Pi = zi(oVq + ln(Li + p) — lni?i H----— (L\ + p — R\)d + 0(d 2 ) 

KB-!- Z2 

is the difference between the boundary electrochemical potentials for (BVi), and 

[»] _ z\{L\ + p — R\) [i] _ z\{L\ + p — R\) z\(L\ + p + R\) 

0 ln(Li + p) — In R\ 1 ln(Z*i + p) — In R\ 2 


Similarly, the flux J-J^ is given by 


j}° ] = 11 + 


2{\z 1 -z 2 ) [o]\ Diw [ ° ] p d i J 
Z1Z2 1 ) Z\H(1) k B T 


+ 0(d 2 ), 


where 


Pi° = ziCqVq + InLi - ln(i?i + p) + — - ^ + (L\ - R\ - p)d + 0(d 2 ) 


k B T' ^ ... Z2 

is the difference between the boundary electrochemical potentials for (BVo), and 


[o] _ zi(Li - Ri - p) [ 0 ] _ zi(Li - Ri - p) _ zi(Li + Ri+ p) 

0 In L i — ln(i?i + p) 1 In L\ — ln(i?i + p) 2 

Recall that s = L\jR\. One gets 


Pi =zi(oV 0 + ln.s + ^(s _ 1 )R\d + 0(pd , d 2 ), 


k B T 


s — 1 

Wn' =w- Z\Rl + 


sRl Z 2 

s In s — (s — 1) 


w\ =- 


In s 
s — 1 
In s 


s(ln s) 


ziRi - S -tt l z 1 Ri + 


Z!p + 0(p ), 

/sins - (s - 1) 1\ _ , 

V 2.) Zlp + °( p 


t^P i’° =z 1 &V 0 + ]ns- |- + 2Z2 (1 + A)Zl (s - l)Ri d + 0(pd,d 2 ), 

k B l Ri Z2 

fol S - 1 „ s - 1 - Ins _ . 2 x 

ziRi H-^— z\p + 0(p ), 


W ° Ins 

[o] _ s ~ 1 
Wl Ins 


(In s) 

z\Ri - ^l-^i + 


s — 1 — In s 1 


(Ins) 


- 2 ) zip + 0(p 2 ). 


Thus, 


=1 + 7 /w li - +0(p 2 ,pd,d 2 ) t 

p d {° s(zi( 0 V 0 + lns) Ri 

2s In s — s 2 + 1 p 2 [i] [o] 2slns-s 2 + l 2 n 

— =H-7-—-—+0(p), = ^ -zip + 0(p). 


w. 


[o] 

0 


s(s — 1) Ins Ri 


s(lns) 
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Therefore, the ratio between the fluxes ,/f^ and is given by 
j'l _Z\Z 2 + 2(\Z! - z 2 )w^d v$ n S { 1 


J 1 °^ ziz 2 + 2(Xzi - z 2 )w^d Pi° 


=1 + 

As a consequence, 
So (p) = 1 + 


ziz 2 

2s In s — s 2 + 1 


+ 0{d 2 


= |1 + S^E) (m l.] + 0{S) 

/ 1 Mi 


+ 


s + 1 


s(s — 1) In s ^(ziCo^b + Ins)/ R\ 


P 


— +0(p , pd, d ). 


2s In s — s 2 + 1 


s + 1 


.(. - i)m. + . (aW + M ) k and Sl<0) = a (7 W) 

Note that So(p) obtained here agrees with that in (7.16), as expected. The conclusion 
S'i(O) = 0 implies that there is not linear term in d in the approximation for 


7.2.3 Approximation of 

The superposition of (7.16) and (7.19) then gives 


./[^ _ P 2s In s — s 2 + 1 s + 1 

j|°] + V s(s - 1) Ins + s^iCo^o + Ins) 

+ O(p 2 ,d 2 ,pdQ 0 ). 


_P_ 

Ri 


+ gi(V 0 ,s,a,(3) 


pQ 0 

R 2 


(7.20) 


Finally, formula (5.14) in Theorem 5.1 can be obtained from the formula (7.2) by 
multiplying the three factors estimated in (7.4), (7.5) and (7.20). 


7.3 Derivation of formula (5.18) 

We now derive formula (5.18) for the case where approximate electroneutrality bound¬ 
ary conditions are assumed; that is, we assume the boundary concentrations L\ and 
R± for the main ion species and those L 2 and R 2 for the counterion species stay the 
same when tracer is added. Therefore, the electroneutrality boundary conditions are 
only approximate with the tracer concentration p being small. In this approximate 
electroneutrality boundary conditions, there will be boundary layers to compensate 
the imperfect electroneutrality. 

More precisely, for (BVi) with tracer concentrations c^’^(0) = = p > 0 at 

x = 0 and cf’^(l) = = 0 at x = 1, one has 

(BVi): zi(Li + p) + z 2 L 2 = zip « 0 and z\R\ + z 2 R 2 = 0, 

and there will be one boundary layer at x = 0. 

For (BVo) with tracer concentrations c^’^(0) = = 0 at x = 0 and c^’°\l) = 

R^'°^ = p > 0 at x = 1, one has 

(BVo): z\L\ + z 2 L 2 = 0 and z±(Ri + p) + z 2 R 2 = z\p ~ 0, 
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and there will be one boundary layer at x = 1. 

For (BVi), let (<Ado, c iio’ c^oo, <4* lo) he the corresponding values of (<jM, c^, cf’4 c^] 
at the limiting point of the boundary layer at x = 0. Then, it follows from [56] (Lemma 
3.2 and Proposition 3.3 in [56]) and (5.17) that 

=^0 - zrz ~-v In ~ Z2l _; 2 = P 0 + -rf + 0(p 2 ), 

Co(zi-Z2) zi(Li + p) Co(^l ^2) L\ 


_W _r I L 1 

h 00 1 ' Ll + p 


z l~ z 2 


= L] - 


z 1 


[t,i] 

C 1,TO =P 


L\ 

L\ + p 
Li 


zi - z 2 

,2 


= P 


^= l ’Kl 1 + p 


= L 2 — 


P + O(p) 2 , 

Zl (r + O(p) 3 

P + 0{p) 2 . 


(7.21) 


zi - z 2 Li 
z\L 2 


(zi - z 2 )Li 


It is easy to check that zi(c[4o + c^) + z 2 c^ 00 = 0, which is known to be true for 
the limiting point ([56]). 

Similarly, for (BVo), let (</>to, c^, cf’^, c^) be the corresponding values of 
{$-°\ c i°^ c i’°\ c 2 °b the limiting point of the boundary layer at x = 1. Then, 
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In ■ 


—z 2 R 2 
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zi - z 2 Ri 
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P + O(p) 2 , 

p ~ +o ( P ) 3 , 

p + 0(p) 2 . 


(7.22) 


(zi - z 2 )Ri 


One can derive the formula (5.18) following exactly the same procedure as in 
Sections 7.2. Another approach is to make use of the derivation in Section 7.2 with 
necessary modifications. We have used both approaches and, not surprisingly, got 
the same approximate formula (5.18). The latter is simpler and is presented below. 
More precisely, we need some modifications in the derivation in Section 7.2. 

For (BVi), the boundary conditions should be replaced by 


r 

'00 > 


> c i, J to> c i,to> c 2 , J oo) at a; = 0; (0, R u 0, R 2 ) at x = 1. 


For (BVo), the boundary conditions should be replaced by 

(Vo, hi, 0, L 2 ) at x = 0; (^, cf^, cj^lo) at x = 1. 

With these modifications, it is not hard to get, from (5.6) that, for the boundary 
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condition associated to (BVi) in (5.17), 


p[\o + -zi,di) — + ^ 2(^00 + *) + (1 + A)c^oo) + ^(^ 2 ) 

= ~iCoTo + - ~Yj —^ d(^2L\ + (1 + A)L 2 ) + 0(d 2 ,dp), 

pf(l;zi,di) =d(4i? 1 + (l + A)i? 2 ) +0(d 2 ); 
and, for the boundary condition associated to (BVo) in (5.17), 

Pi ^(0; zi, d±) =ziCo^o + d{2L\ + (1 + A)L 2 ) + 0(d 2 ), 

Pi\ 1 i z u di) =zi(o</>£! + ^(2(4,L + 4,4) + (1 + A)4,L) + 0{d 2 ) 

=---——f- d(2R\ + (1 + A)i? 2 ) + 0(d~, dp). 

zi — z 2 /t-l 

Thus, it follows from the definitions of the factors in (7.2) that 
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1,00 

x ( 1 + (zl-z'li; + ( 2(Ll -^ + (1 + A)(L2 - R ^ d ) 

+ 0{d 2 ,dp, p 2 ) 
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Also, Fq in (7.15) should be replaced by 
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A straightforward calculations gives 


(7.23) 


(7.24) 
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where /i(s) is defined in (5.8). 

The most complicate term is Fq/F\ in (7.15). By a close examination of the terms 
involved in the expression, the following modifications are needed. One is to replace 
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p with —Z 2 p/(z\ — Z 2 ). The other is, in the last expression of Fq/F, in (7.15), to 
replace Vo by ^ in the first term and to replace Vo by Vo — 4$ in the second term. 
The combined effect of the latter is the extra term 


{E)= _ ^(1-jjH) _ *2 P 

(*1 - *2) (Li - J^r 2 p - Ri) (ln(L! - ^p) - In R 1 ) *1 ~ ** sR 1 

z 2 P 
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Using the expressions of A$, Aq 0 ^, and A^A?^ in (7.17), one gets 


(E) = g 2 (s,a,P)- 


-Z2 p 
’ Z\ - Z2 R‘\ ' 

where g 2 (s,a, (3) is given in (5.9). Therefore, we have 
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and hence, 
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+ (gi(Vo,s,a,P) + g 2 (s, a, (3)j — ^ 2 ° + 0(p 2 , d 2 , pdQ 0 ). 


(7.25) 


The formula (5.18) in Theorem 5.2 then follows from (7.2), and the estimates 
(7.23), (7.24) and (7.25). 
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